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RESUMEN 


El Metodo del Elemento Finito (MEF) ha llegado a ser una herramienta muy poderosa en la 
solucion numerica de un amplio numero de problemas en ingenierla. Las aplicaciones van 
desde el analisis de esfuerzos y defonnaciones en aeronaves, edificios, barcos, automoviles, 
puentes, ademas de problemas de flujo de calor en cuerpos, de fluidos, campos magneticos, 
filtraciones, etc. Con el avance de la tecnologla de los procesadores pueden modelarse 
problemas complejos con relativa facilidad. En la actuabdad existe un gran numero de 
programas comerciales y de investigacion cuyo planteamiento es con base en el MEF. 

En este Metodo una region compleja que define un medio continuo se discretiza en formas 
geometricas simples llamadas elementos finitos. Sobre esos elementos se consideran las 
propiedades del material y las relaciones gobemantes, expresadas en terminos de valores 
desconocidos en las fronteras del elemento. El sistema de ecuaciones de equilibrio se logra 
cuando se ensamblan las cargas y condiciones de apoyo. La solucion de esas ecuaciones de 
equilibrio dan lugar a una solucion aproximada en el continuo. 

OBJETIVO GENERAL 

Establecer en fonna clara y completamente didactica los fundamentos matematicos del 
metodo del elemento finito con base en la teoria elastica lineal para que la poblacion 
profesional de ingenieria civil aplique correctamente por medio de programas comerciales 
la teoria. 

OBJETIVOS ESPECIFICOS 

- Desarrollar la teoria del Elemento Finito mediante la fonnulacion integral al resolver 
numericamente ecuaciones diferenciales. 

- Desarrollar tecnicas numericas utilizando los sistemas coordenados para la evaluacion de 
integrales cuya solucion analitica es imposible. 

- Resolver por medio del elemento finito los diferentes elementos unidimensionales y 
bidimensionales. 
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CAPITULO I 

INTRODUCTION 

1.1 ASPECTOS HISTORICOS 

Las ideas basicas del Metodo se originaron gracias al avance en el analisis de estructuras de 
aeronaves. En 1941, Hrenikoff presento una solucion de problemas de la elasticidad usando 
el “metodo de trabajo del marco”[l]. En un artlculo publicado en 1943 [2], Courant uso 
interpolacion polinomial por partes sobre subregiones triangulares para modelar problemas 
de torsion. Turner y otros investigadores obtuvieron matrices de rigidez para armaduras, 
vigas y otros elementos y presentaron sus hallazgos en 1956. Clough, en 1960, fue el 
primero en emplear el tennino de elemento finito[3] . 

En los primeros anos de la decada de 1960 los ingenieros usaron el metodo para obtener 
soluciones aproximadas a problemas de analisis de esfuerzos, flujo de fluidos, transferencia 
de calor y otras areas. En 1955 Argyris cimento el analisis del elemento, publicando los 
teoremas de energia y metodos matriciales[4]. El primer libro sobre elementos finitos fue 
publicado por Zienkiewicz y Cheng, en 1967[5], A fines de la decada de 1960 y principios 
de la siguiente, el analisis por elemento finito se aplico a problemas no lineales y de 
grandes deformaciones. Oden publico su libro sobre aplicaciones no lineales en 1972. Las 
bases matematicas se establecieron en la decada de 1970. 

Actualmente se cuenta con una gran cantidad de programas comerciales de elemento finito, 
tales como: SAP 2000, Ansis, Adina, etc. 


[1] A.Hrenikoff, ‘Solution of problems in elasticity by the framework method’, J. Appl..Mech.,A8, 169-75, 1941. 

[2] R. Courant, ‘Variational methods for the solution of problems of equilibrium and vibration’, Bull.Am.Math.Soc,49,l- 
23,1943. 

[3] R.W. Clough, ‘The finite element in plane stress analysis’, Proc.2nd A.S.C.E. Conf. on Electronic Computation, 
Pittsburg, Pa., Sept. 1960. 

[4] J.H. Argyris, Energy theorems and structural analisis, Butterworth, 1960 (reproducido de Aircraft Eng., 1954-55). 

[5] O.C Zienkiewicz, Y.K. Cheng ‘The finite element method, (1967). 
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1.2 SOLUCION DE PROBLEMAS CON VALORES EN LAS FRONTERAS 

Existen varios metodos numericos para resolver ecuaciones diferenciales: 

1) Metodo de Diferencias Finitas (MDF) 

2) Metodos Variacionales 

3) Metodos Integrales o Elementos de Frontera (MEF) 

4) Metodo de Elemento Finito (MEF) 

4.1) Metodo basado en una formulacion integral (Residuos Pesados) 

4.2) Metodo basado en una formulacion energetica (Formulacion General) 

1.2.1 METODO DE DIFERENCIAS FINITAS 

Este metodo aproxima las derivadas en la ecuacion diferencial usando ecuaciones en 
diferencia o en incrementos. Se utiliza para resolver problemas de transferencia de calor, 
mecanicos, de fluidos y problemas bidimensionales con fronteras paralelas a los ejes 
coordenados. El metodo es dificil de aplicar a una region curva o fronteras irregulares. 

1.2.2 METODO VARIACIONAL 

Involucra la integral de una funcion que produce un numero, cada nueva funcion produce 
un numero mas bajo. La funcion que produce un numero mas bajo, tiene la propiedad 
adicional de satisfacer una ecuacion diferencial especifica. El Metodo Variacional, es la 
base de muchas formulaciones del elemento finito, pero tiene la desventaja, que no es 
aplicable a ninguna ecuacion diferencial que contenga un tennino en la primera derivada. 
En ingenieria civil las ecuaciones de elementos mecanicos estan gobemadas por ecuaciones 
diferenciales de segundo grado y en el caso de placas de cuarto grado. 
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1.2.3 METODO DE ELEMENTO FINITO (MEF) 

Es un procedimiento numerico para resolver problemas flsicos gobemados por una 
ecuacion diferencial o un teorema de energla. Dos de las caracterlsticas que lo distinguen de 
otros metodos son: 

a) El metodo utiliza una formulacion integral para generar un sistema de ecuaciones 
algebraicas. 

b) El metodo utiliza funciones continuas para la aproximacion de una cantidad o 
cantidades desconocidas. 

La segunda caracteristica distingue al MEF de otros procedimientos numericos que utilizan 
una formulacion integral. El MEF se puede dividir en cinco partes: 

1) Se discretiza la region: Esto incluye la localizacion y numeracion de puntos 
nodales, asi como especificar el sistema coordenado (topologia Fig. 1.1). 




MAYOR ANCHO OPTIMO 

DE BANDA 

Fig. 1.1 Topologia de un modelo marco. 

2) Especificar la ecuacion de aproximacion 
Se debe especificar el orden de aproximacion lineal o cuadratica y la ecuacion debe de 
escribirse en tenninos de valores nodales desconocidos. Se escribe una ecuacion para cada 
elemento. 
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3) Desarrollar el sistema de ecuaciones usando el metodo de fomiulacion integral 
de Galerkin (Residuos Pesados) o la transformacion energetica. 

4) Resolver el sistema de ecuaciones. 

5) Calcular las cantidades de interes o incognitas, que usualmente son las derivadas 
de los parametros que incluyen las componentes de esfuerzos, defonnaciones, 
etc. 

Existen dos enfoques para abordar la solucion de problemas en el MEF : 

a) Fomiulacion Integral 

b) Fomiulacion Energetica 

1.3 FORMULACION INTEGRAL, METODO DE RESIDUOS PESADOS 

Los metodos numericos mas usados para resolver ecuaciones diferenciales ordinarias son 
los conocidos como Metodo de Residuos Pesados. Estos metodos empiezan por suponer 
una solucion a la ecuacion diferencial que satisfaga las condiciones de frontera. Dicha 
solucion contendra ciertos parametros de los cuales pueden ser ajustados para minimizar los 
errores. 

Existen diferentes metodos para minimizar el error o residuo, tal que dicha solucion 
propuesta se aproxime tanto como sea posible a la solucion exacta. Los Metodos de 
Residuos Pesados son: 

1) Metodo de Colocacion Puntual. 

2) Metodo de Galerkin. 

3) Metodo de Minimos Cuadrados. 

4) Metodo del Subdominio. 


[Referencia] ‘Analisis Avanzado de Estructuras’, apuntes del Ing. Hugo Hernandez Barrios. 

1.3.1 PROCEDIMIENTO 

i) Se tiene una ecuacion diferencial. 

&£> = F en Q de la ecuacion diferencial (1.1) 
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donde: 


8 es el operador diferencial 


d 2 d_ 
dx 2 dt ’ 


O es una variable \u,v,w, 

F es una funcion. 

Cl es el dominio de la region 


ii) Se supone una solucion a la ecuacion. 

0*«0 ( 1 . 2 ) 

donde 0 es una aproximacion a la solucion cerrada, 0 , tal que se puede escribir como un 
polinomio, 

®'=t. c Xi <>.3) 

(=1 

0* = f 0 (x) + CJ\ (x) + C 2 f 2 (x) + . . . (1.4) 

donde / (x) puede ser cero. 

La solucion supone que las funciones, satisfacen las condiciones de frontera de 

la ecuacion diferencial y contiene un nurnero “indetenninado” de parametros C 1 ,C 2 .... 


iii) Se sustituye 0 en la ecuacion diferencial, 


£0* =F 


en Cl 


(50* -F = 0 

solo si 0* es igual a la solucion exacta, en otro caso: 


(1-5) 

(1.6) 


<50* -F = R (1.7) 

donde R es el Residuo o error 


iv) Se encuentran los valores de C, , C 2 . . . incognitas que minimi z an el error ( R ) en el 
dominio de la solucion A<x<B (condiciones de frontera). 
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v) Se pesa o pondera el error en todo el dominio de la funcion por medio de una 
funcion de peso, Wj. 

J[^D’-F]^£/Q = 0 (1.8) 

Q 

es decir: 

fRW i dQ = 0 (1.9) 

Q 

donde el error o residuo es: 

/? = [<5S>*-f] (1.10) 

La mejor aproximacion de la solucion de la ecuacion diferencial, minimiza el error, R. 
Dependiendo del peso W; que se proponga, sera el metodo utilizado, asl, 

a) Se conoce como Metodo de Colocacion Puntual si la funcion de peso es: 

W t =S(x t ) (1.11) 

donde 5 (x i ) es la funcion Delta de Dirac evaluada en un punto en el dominio. 

Por definicion: 

*2 

J S^x^G^x^dx = G(x ; ) en donde: x 1 < x j < x 2 (1.12) 


b) Se conoce como el Metodo de Galerkin si la funcion de peso es: 

W i =f i (1.13) 

donde f] es un polinomio (una funcion) del mismo orden que la funcion de aproximacion 

o . 


c) Se conoce como el Metodo de Minimos Cuadrados si, 

W i =s = R 


la funcion de peso es el mismo error, el cual se minimiza: 


d 

dC \ 


js-sdQ 


0 


(1.14) 


(1.15) 
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haciendo 


se tiene un nuevo error. 


se minimi z a el nuevo error. 


5 


dC, 


J[r] 2 JQ 


0 


Q. 


ds- R 
~dC~ 


0 


( 1 . 16 ) 


( 1 . 17 ) 


( 1 . 18 ) 


d) Se conoce como el Metodo de Momentos o del Subdominio si, 

W t = x J 


la funcion de peso es un polinomio: x J . 


( 1 . 19 ) 


EJEMPLO 1.1 

Resolver la siguiente ecuacion diferencial: 

dx 

— = -x 
dt 

Usando: 

a. El metodo de colocacion puntual. 

b. El metodo de Galerkin 

c. El metodo de momentos o del subdominio 
La condicion inicial es: 

x = 1 en t = 0 


en 0 < t < 1 


I. Se propone como solucion aproximada un polinomio de segundo grado, 

x* = 1 + Cj • t + C 2 • t 2 

que satisface: 

x* = 1 en t = 0 

II. Se sustituye en la ecuacion diferencial original, 

dx * 


dt 


- + x =R 


( 1 . 20 ) 


( 1 . 21 ) 

( 1 . 22 ) 

( 1 - 23 ) 


( 1 . 24 ) 
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^[l + C, ■t + C 2 -t 2 ~\ + [l + C l -t + C 2 -t 2 ] = R (1.25) 

desarrollando, 

C x +2C 2 -t + \ + C x -t + C 2 -t 2 =R (1.26) 

tomando terminos comunes, el error es: 

C 1 [l + t] + C 2 [2t + t 2 ] + l = R (1.27) 

III. Se minimiza el error: 

J RW i dQ = 0 (1.28) 

n 

a) METODO DE COLOCACION PUNTUAL 

El dominio de la ecuacion diferencial se muestra en la figura 1 .2: 


I 1 1 1 

0 0.3 0.6 1 

Fig. 1.2 Dominio de la ecuacion diferencial del ejemplo 1.1. 


Existen dos incognitas o parametros desconocidos: C/ y C 2 por lo que se necesitan dos 
valores de peso, 


W l =£(0.30) 

W 2 =S( 0.60) 

Evaluando las ecuaciones (1.29) y (1.30) en el dominio de la ecuacion diferencial, 

1 

J[C, (l + t) + C 2 (2t + t 2 ) + l] -^(0.30) = 0 


1 

J[c, (l + f) + C 2 (2t + r) + l] -^(0.60) = 0 


(1.29) 

(1.30) 


(1.31) 

(1.32) 


por definicion: 

*2 

jS(x i )G(x)dx = G(x i ) (1.33) 

xl 


tal que la solucion de las integrales son: 
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Cj (l + 0.30) + C 2 ^2(0.30) + (0.30) 2 j + 1 = 0 

(1-34) 

C, (1 + 0.60) + C 2 (2 (0.60) + (O.6O) 2 ) + 1 = 0 

(1.35) 

1.30Cj + 0.69C 2 +1 = 0 

(1.36) 

1.60Cj +1.56C 2 +1 = 0 

(1-37) 

Resolviendo el sistema: 


C, =-0.9414 

(1.38) 

C 2 = 0.3245 

(1-39) 


Sustituyendo en la solucion propuesta: x* = 1 + C x ■ t + C 2 ■ t 2 , la solucion aproximada es: 

x* = 1-0.9414? + 0.3245? 2 (1.40) 


b) METODO DE GALERK1N 


Se necesitan dos funciones de peso ya que existen dos parametros incognitas, Cj y C 2 , se 
proponen: 

W x =t (1.41) 

W 2 =t 2 (1.42) 

En donde al menos un polinomio es del mismo orden que la funcion de aproximacion, en 
este caso de segundo orden. 


J R -Wd Q = 0 

Q 


J[c,(l + t) + C 2 (2t + t 2 ) + l 


•tdt = 0 


0 

1 

J[c,(l + t) + C 2 (2t + t 2 ) 

0 


+1 


■t 2 dt 


0 


Resolviendo las integrates: 


J[c 1 (t + r) + c 2 ( 2 t 2 


+ t 3 ) + t 


ut = 0 


(1.43) 

(1.44) 

(1.45) 


(1.46) 
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c, 


c, 


t 2 t 3 


- 2t 2 

1 

+ C 2 

— + ■ 

2 3 _ 


_ 3 

"i r 


'2 r 

— 1 — 

+ c. 

— 1 — 

_2 3_ 

z 

_3 4_ 


r 

+ — 
2 


+ - = 0 
2 


= 0 


0.8333C, +0.9166C 2 +0.50 = 0 


Para la segunda integral: 


|[Cj [t 2 +t 3 ) + C 2 (2t 3 +t 4 ) + t 2 


Q 

c, 


V 3 

t 4 1 

+ C 2 

- 2t 4 

t 5 

t 3 

— + — 


+ — 

+ — 

3 

4 


4 

5 

3 


|<7t = 0 


= 0 


1 

ll 

+ C 2 

n 

ll 

- + — 

-+- 

3 

4 1 


12 

5J 


+ - = 0 
3 


0.5833C,+0.70C 2 +0.3333 = 0 
Resolviendo simultaneamente (1.49) y (1.53): 


(1.47) 

(1.48) 

(1.49) 

(1.50) 

(1.51) 

(1.52) 

(1.53) 


0.8333^+0.9166^+0.50 = 0 (1.49) 

0.5833C, + 0.70C 2 + 0.3333 = 0 (1.53) 

Se tiene, 

C, =-0.9145 (1.54) 

C 2 = 0.2859 (1.55) 

Sustituyendo en la solucion propuesta: x* = 1 + C x ■ t + C 2 • t 2 la solucion aproximada es: 


x* = 1 - 0.91456t + 0.2859t 2 (1 .56) 

c) METODO DE MOMENTOS O SUBDOMINIO 
La fimcion de peso es: W f = x j , se proponen: 
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w { =r 
W 2 =t 3 

por existir dos parametros C; y C 2 , incognitas. Por lo que en: 

J R W t d Q = 0 

Q 

1 

J[c 1 (l + t) + C 2 (2t + r) + l’ 

0 

1 

J[c 1 (l + t) + C 2 (2t + r) + l’ 

0 

J[c,(t 2 +t 3 ) + C 2 (2t 3 +t 4 ) 


Para la primera integral: 


■Vdt = 0 


■rdt = 0 


+ r 


\dt = 0 


C, 


ft 3 

t 4 l 

+c, 

IV 

t 5 l 

t 3 

-+— 

1 

+— 

L3 

4 

z 

2 

5 _ 

3 


= 0 


realizando operaciones: 


C, 


n 

ll 

+c 2 

n 

ll 


-+- 

L 3 

4 J 


L2 

5 J 


+ - = 0 
3 


0.5833C, +0.70C +0.3333 = 0 


Para la segunda integral: 


1 

J [Cj (\ + 1) + C 2 (2t + 1 2 ) + l]-t 3 dt = 0 


}[c i (t 3 +t 4 ) + C 2 (2t 4 +t 5 ) 


+ t 


\dt = 0 


C, 


ft 4 

t 5 l 

+ c 2 

r 2t 5 

t 6 i 

t 4 


1 — 

H 

L4 

5 


5 

6 

4 


= 0 


(1.57) 

(1.58) 

(1.59) 

(1.60) 

(1.61) 

(1.62) 

(1.63) 

(1.64) 

(1.65) 

( 1 . 66 ) 

(1-67) 

( 1 . 68 ) 
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n 

ll 

+c. 

r 2 

ll 

-+- 

-+- 

L4 

5J 

z 

L5 

6 J 


realizando operaciones: 

0.45C,+0.5666C 2 +0.25 = 0 


Resolviendo simultaneamente: 


(1.69) 


(1-70) 


Se tiene: 


0.5833^+0.70^+0.3333 = 0 

(1.71) 

0.45Q+0.5666C 2 + 0.25 = 0 

(1.72) 

Cj =-0.8935 

(1-73) 

C 2 = 0.2684 

(1.74) 


x* = 1 - 0.8935? + 0.2684? 2 


RESUMEN: 

La ecuacion diferencial: 


acepta la solucion: 


Tal que: 


Sustituyendo en la ecuacion (1.76), 


dx 

hX = 

dt 


0 


x = e 


-t 


dx 

— = -e 
dt 

-e~* +e~‘ =0 


Por lo que la solucion exacta de la ecuacion diferencial (1.76) es: 


x = e 


Las soluciones numericas anteriormente calculadas son: 


Metodo de Colocacion Puntual: 
Metodo de Galerkin: 

Metodo del Subdominio: 


x = 1-0.9414? + 0.3245L 


x = l-0.91456? + 0.2859r 


x = 1 - 0.8935? + 0.2684?* 


(1.75) 

(1-76) 

(1-77) 

(1-78) 

(1-79) 

(1.80) 
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En la Figura 1.3 se observan las graficas de las soluciones numericas en el intervalo de 
solucion 0 < t < 1 . 



0.5 

t 


Fig. 1.3 Comparacion de resultados de los metodos de residuos pesados. 


COMENTARIOS: 

I. Para este problema en particular, el metodo de colocacion puntual y el de Galerkin se 
aproximan mejor a la solucion exacta que el del Subdominio. 

IF El metodo de colocacion puntual, depende de los puntos que se elijan para la 
solucion y que ademas deben pertenecer al dominio de la ecuacion. 

III. El metodo del Subdominio tiene la desventaja que no se tiene control, en caso de 
programarse, del orden del polinomio. 

IV. El metodo de Galerkin, es el mas facil de programar ya que se puede colocar la 
funcion de peso igual al polinomio de aproximacion. 
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EJEMPLO 1.2 

Resolver la siguiente viga usando: 

a) Metodo de colocacion puntual 

b) Metodo de galerkin 

c) Metodo de mlnimos cuadrados 

La ecuacion que gobiema el momenta flexionate en una viga es: 

(LSI) 

Donde y es la flecha o defonnada de una viga. Para el caso en particular de la viga 
simplemente apoyada sometida a un momenta flexionante, M o , en sus extremos, como se 
observa en la Figura 1.4, las condiciones de frontera son: 

En x = 0 : y - 0 (1-82) 

En x = L : y = 0 (1.83) 


Mo 


Mo 


(7 


Mo 


7 


El 


E> 


Mo 


Diagram a de Moment os 

Fig. 1.4 Diagrama de momentos flexionantes del ejemplo 1.2. 


La solucion exacta del problema es la ecuacion: 

Mox , y- \ 

' = wF ( x - l) 

De la ecuacion diferencial (1.81) se tiene que la ecuacion a resolver es: 


(1.84) 


(1.85) 


Se propone una aproximacion a la solucion del tipo: 
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Que satisface las condiciones: 


* . nx 

y , , = Asen — 

* 0) L 

(1.86) 

o 

II 

o 

II 

H 

(1.87a) 

x-L: y = 0 

(1.87b) 


Sustituyendo M (x , = M o en la ecuacion (1 .85), la ecuacion diferencial a resolver es: 


Mo-EI^- = 0 


dx 


Sustituyendo la aproximacion a la solucion: 


Mo-EI d -^ = R 
dx 


Donde segun (1.86) 


. . nx 

y {x) = Asen — 

es la aproximacion a la solucion. Sustituyendo en (1.89), 

d 2 


Mo - El - 


dx~ 


. nx 
Asen — 
L 


= R 


Es necesario calcular las primeras dos derivadas de (1.86): 

dy( x ) , n nx 

dx L L 



dx 1 


. n nx 
-A -—sen — 


V 


L 


Sustituyendo en (1.90) se tiene: 


Mo + EI 


An 2 

L 2 


nx 
sen — 
L 


R 


(1.88) 

(1.89) 

( 1 . 86 ) 

(1.90) 

(1.91) 

(1.92) 

(1.93) 


a) METODO DE COLOCACION PUNTUAL 


Se pesa el error, 

jRW ; d Q = 0 

Q 

el dominio de la ecuacion es: 
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H 


I 

0 L 

Fig. 1.5 Dominio de la ecuacion del ejemplo 1.2. 


La funcion delta de Dirac se evaluara en x = ~- Se tiene un solo parametro incognita “A”. 
Sustituyendo: 


resolviendo la integral: 


Despejando: 


L 

r 

An 2 nx 



[ 

Mo + El — — sen — 

5 


J 

0 

1 

pi 

1 


uJ 


dx = 0 


Mo + El — — sen — 
L 2 L 




V ^ / 


= 0 


Mo + El lEL- sen — = () 
L 2 2 

An 1 

Mo + El — — = 0 
IS 


A = - 


MoL 1 


k 2 EI 


Sustituyendo en la solucion aproximada, (1.86): 


F ( 


*) 


MoL 


t: 2 EI 


2 A 

sen 


m 

L 


(1.5) 


b) METODO DE GALERK1N 

Se pesa el error, 


J R W t d , Q = 0 


La funcion de peso, W i , es de la misma fonna matematica que la funcion de aproximacion, 

jjr nx 

W ; = sen — 


Sustituyendo: 
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0 L 


,, nr An nx 

Mo + El — 7 — sen — 
L 2 L 


sen—dx = 0 
L 



sen — dx + El 
L 


An 


2 L 


L 1 


sen- 


nx 

L 


dx = 0 


Para resolver las integrates, se recurre a las siguientes formulas: 


r , , cos (/roc) 

sen ( mil dx = - 

J m 

f sen 2 u ■ dx = ~ — sen2u + — 

J 4 2 


si: 


nx 
u = — 


L 


du = — dx 
L 



—du = dx 
n 


Mo 


L nx 

cos — 

n L 


L 


+ EIA 

o 


n^L 

L 2 n 


nx 

1 2 nx t 

— sen 

4 L 2 


L 

0 


0 


MoL 


-cos 


n 


r nx'' 

L , n 


+ EIA — 


0 L _ 


1 2nx nx 

- — sen 1 

4 L 2L 


0 


-Mo — [cos (n)- cos (0)] + £Z4 


n 


-—sen^-(L) +—sen(2n(0)) + 

4 L y ’ 4 v 2 L 


(L) n( 0) 


2 L 


0 


-Mo-[-l-l] + EIAj 


— sen2n +— sen 

4 4 


( 0 ) + 


n 

~2 


0 


2 Mo- + EIA- 
n L 


n 


= 0 


Despejando la incognita: 


A 


-2 Mo- 
n 


2 L 
Eln 2 
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Se tiene, 


AMo 1} 

~eT 1? 


Sustituyendo en la solucion aproximada (1.86): 




4 Mo L 2 
~EJ^? 


nx 
sen — 
L 


(1.6) 


c) METODO DE MINIMOS CUADRADOS 


Se trata de encontrar un nuevo error: 


ei? = J[i?(x )] 2 


Asl, el nuevo error es: 


eR = \ 


-• r nr j KX 

Mo + El A — — sen — 


V 


L 


dx 


Desarrollando: 


L 

_ f 

9 n nx 

n 2 

2 

2 nx 

ei? = 

Mo~ + 2 Mo EIA —r sen 1- 

EIA — 

sen — 

J 

0 

1 

L 2 

L 


\mo 2 x + 2MoEIaE 


nx 

-cos — 

+ 

KTA 7 * 2 

EIA — 

2 

~L~ 

1 2 nx nx 

— sen 1 

) 

L 


L 


L 


_n _ 

4 L 2L 

J 


R = < J Mo 2 -2MoEIA— [-cos;r-cos(0)] + 
L 


EIA 


n 

~L 2 


- — sen2n + — sen (0) + 
4 4 w 


9 n 

2 1 2 

n 

Tz" 

n 

mo 2 l+amoeia~+ 

EIA — 



L 

L 2 _ 

_;r_ 

_~ 2 _ 


e R = Mo 2 L + AMoEIA — + 

, r 2 I 2 

\l 

EIA 

L 

l 2 _ 

_2_ 


El nuevo error es: 


g R = Mo 2 L + A [4 MoEIn] + A 2 


~L~ 

n 2 


El — 

_n _ 

L 2 _ 
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Minimizando el nuevo error con respecto a “A 


d eR 
dA 


0 


d g R 

- JLJ Mn 2 I 4 - A 

AMoEI/r 

+ A 2 

~ L 

' Eln 2 1] 

dA 

— A 1V1U J-J 1 /I 

5^1 

L 


i 2 J J 


AMoEIn 

L 


+ AL 


Eln 2 
L 2 


= 0 


= 0 


despejando, 


A = -AMo 


L 2 


Eln 3 


y sustituyendo en la solucion aproximada ( 1 . 86 ): 

* 4Mo Is nx 




El n' 


-sen- 


L 


(1.7) 


RESUMEN: 


Resumiendo, las soluciones son: 




Solucion exacta: 


Mox / r . 

= (x L) 

2EI V ’ 

(1.84) 



Mo Is nx 


Metodo de Colocacion Puntual: 


= r — sen — 

n-EI L 

(1.5) 

Metodo de Galerkin: 


AMo Is nx 

( 1 . 6 ) 

yU 

= sen — 

El n l L 

Metodo de Minimos Cuadrados: 

* 

AMo Is nx 

= 5 - sen — 

El n 3 L 

(1.7) 


Comentarios: 

I. Para este caso, no es posible establecer cual metodo es el mas exacto, el error 
depende de la funcion de aproximacion y de la ecuacion resuelta. 

II. Los puntos de colocacion o subregiones seleccionadas afectan la exactitud de la 
solucion. 

III. Lo importante de esta seccion es ver que la solucion numerica de ecuaciones 
diferenciales se pueden fonnular en terminos de una integral. La formulacion integral 
es una caracteristica del MEF. 
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CAPITULO II 

CONCEPTOS BASICOS 

2.1 ELEMENTO LINEA 

El elemento linea (metodo de interpolacion) es un segmento de llnea con longitud L y dos 
nodos, uno en cada extremo. Los nodos son denotados por i y j y los valores nodales por 
® ; y ® . . El origen del sistema coordenado esta en la izquierda del nodo i . 

El parametro® varla linealmente entre los nodos (Fig. 2.1), por lo que la funcion de 
variacion es: 

</> = a l + a 2 x (2-1) 



Fig. 2.1 Variacion lineal entre dos nodos. 


Los coeficientes a, y a 2 se determinan con las condiciones en los nodos. 


Sustituyendo en (2.1), 


(j) = ® ; en x = x ; (2.2) 

^ = ® ; en x = Xj (2.3) 

®,. =aj+a 2 x ; . (2.4) 

® ; . =a l +a 2 x j (2-5) 
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Multiplicando la primera por (-1): 

-0,. = -a x - a 2 x i 
+0 . = +a x + a 2 Xj 

y sumando termino a termino: 

a 2 [*;-*,] = ft;-0, 

Despejando el valor de a 2 : 


De la segunda ecuacion (2.6): 

Despejando a x \ 

Sustituyendo a 2 obtenido en (2.8), 


-ft,- 

Xj ~ x t 


0 v = a, + a 2 Xj 


flj = O . - a 2 x . 


6 bien. 


que tambien se puede escribir: 


Realizando operaciones: 


a { = O. - 


0 -0 

J l 

Xj-x, 


0.A.+0.A. 

a t = 0 1 J -^ — L 


<h =ft ; +- 


Xj-Xi 


ft,*; - ft;*; 


*;“*, 


0 (.V -.V ) • 0 r V,-0 




*; — */ 


Por tanto, 


0 A- - 0 ..V,.- 0 ..V / - 0 ..V > . 

*;“*, 


= 


ft/*; ~ft;*, 

x j -x i 


(2.6) 

(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 

(2-14) 

(2.15) 

(2.16) 
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Sustituyendo el valor de a 2 obtenido en (2.8) y el de a, obtenido en (2.16), en la funcion de 
interpolacion, 

</) = a l +a 2 x (2.17) 


Se tiene, 


De la figura 2.1, 


0 -O 

J 

Xj-x, 

L=Xj- X; 


0,x, - O ,.x 
<j) = ^^ ^ + 


Xj-X, 


La ecuacion (2.18) se puede escribir, 


(2.18) 


(2.19) 


Organizando tenninos, 


$ 


0,JC. -0 ,.JC. +0.JC-O.JC 

‘ J J 1 J * 

L 


</> = &; 



+ 0 

J 


x-x t 

L 


Llamando: 


( 2 . 20 ) 


( 2 . 21 ) 


*■= ’l 

(2.22) 

3: 

11 

X 

t~< 1 

_x 

(2.23) 


donde N i ,N j son llamadas funciones de fomia o de interpolacion. Sustituyendo en (2.21), 


que se puede escribir corno: 


</> = A60.+AL0, 


(h = \N. TV." 

T L 1 j J 


0 

l 

0 


Matricialmente: 


^ = [2V] {0} 


(2.24) 


(2.25) 


(2.26) 
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donde [iV] es la matriz que contiene a las funciones de forma y {®J es el vector que 
contiene los valores nodales de la funcion y (j) es la funcion que interpola los valores 
nodales. 


[Referencia] ‘Analisis Avanzado de Estructuras’, apuntes del Ing. Hugo Hernandez Barrios. 


2.2 PROPIEDADES DE LAS FUNCIONES DE FORMA 

1. Cada funcion de forma tiene un valor de uno en su propio nodo y es cero en 
cualquier otro nodo (Fig. 2.2). 

2. La sunra de las funciones de forma en cualquier punto del donrinio es uno. 

3. Las funciones de forma son siempre polinonrios del nrismo orden que la ecuacion 
de interpolacion original. 

4. La derivada de las funciones de forma respecto a x es cero (en este caso). 


N; = 


Xj -x 


L 



A T X-Xi 

N , = L 

7 L 



Fig. 2.2 Grafica de las funciones de forma. 
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EJEMPLO 2.1 

Se usa un elemento unidimensional para aproximar la distribucion de temperature en una 
estructura. La solucion indica que las temperatures nodales son: 

0,=1OO°C (2.27) 

O j = 85° C (2.28) 

Detenninar la temperatura en un punto a 4 cm del origen y la temperature gradiente en el 
elemento. Los nodos i y j estan localizados a 2.5 y 7 cm del origen, respectivamente (Fig. 
2.3). 



Fig. 2.3 Distribucion de temperatura del ejemplo 2.1. 


Se tiene: 


Por lo que: 


x i = 2.5 cm , Xj = 1cm y x = 4 cm 
L = 7 - 2. 5 = 4. 5 cm 
0, =100 °C 


0 ; =85°C 


Sustituyendo en la ecuacion (2.21): 


<t> 

4 



0 + 


x-x t 

L 


0 


7-4 

4.5 


100 + 


4-2.5 

4.5 


85 


(2.29) 

(2.30) 


(j) = 95° C 


(2.31) 
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El gradiente de temperatura es: 


d(/) _ Q>j 
dx L 


85-100 

4.5 


r 

-3.33° — 
cm 


(2.32) 


2.3 SISTEMAS COORDENADOS 

Las soluciones por Elemento Finito requieren de la evaluacion de integrales. Algunas de 
ellas son faciles de evaluar pero en algunos casos son muy dificiles. Para algunas integrales 
es imposible su solucion analitica tal que, para realizar su evaluacion es necesario usar 
alguna tecnica numerica. 

Las dificultades asociadas en la evaluacion de integrales puede disminuir cambiando las 
variables de integracion, esto involucra escribir la integral en un nuevo sistema coordenado. 
Los sistemas coordenados se clasifican en: 

a) Locales s, q 

b) Naturales ij 

c) De area L 


2.4 SISTEMAS DE COORDENADAS LOCALES 5 , q 


Las funciones de forma lineales obtenidas previamente en la seccion 2.1, estan referidas a 
un sistema de coordenadas ubicado a la izquierda del nodo i (Fig. 2.4). 


x 


xi 


i 


j 



N; 





x-x f 

L 


Fig. 2.4 Funciones de fonna del elemento linea, referidas a la izquierda del nodo i . 
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La desventaja de estas funciones de forma se presenta cuando se tratan de evaluar integrales 


que involucran productos de las funciones de forma, tal como: 


\N i (x)-N j {x)dx 


(2.33) 


| N: (x) dx 


(2.34) 


La integration de estas ecuaciones se facilita desarrollando nuevas funciones de forma 
referidas a un sistema coordenado relativo, cuyo origen esta localizado en el elemento, este 
sistema coordenado es conocido como Sistema Local. Los sistemas locales mas comunes 
para un elemento unidimensional son los que tienen el origen del sistema coordenado en el 
nodo / ; y el que tiene el origen del sistema coordenado en el centro del elemento. 

Las funciones de forma para un sistema coordenado local ubicado en el nodo i (Fig. 2.5), 
se obtiene reemplazando en las funciones de fonna (2.22) y (2.23), 


x = x, + S 


(2.35) 


x 


S 


J 


L 


Fig. 2.5 Sistema coordenado local ubicado en el nodo i . 


Asi se tiene, 



Xj -x 


L 


Xj -Xj-S _ Xj - x . . S' _ | S 

L L L L 


(2.36) 



X-X,. _ x, + S- x, _ s 


(2.37) 


L L L 


Las funciones de fonna obtenidas con respecto al sistema local, S, son: 


iV,(S) = l-| 


(2.38) 
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7V.(S) = | (2.39) 

Que se pueden expresar en funcion de x , 

N f (x) = l-j (2.40) 

N j (x) = j (2.41) 

Estas funciones de fomia son iguales en su propio nodo y son cero en el otro. Su suma es 
uno en todo punto. 


EJEMPLO 2.2 

x i 

Evaluar la integral: J Nf(x)dx 


Para obtener la solucion se cambian los limites de integracion, lo que significa cambiar el 
origen de referencia de la funcion de forma N t (x) — » Af (£) , 


Evaluando la integral: 


j 

| N 2 ( x)dx = 



ds 


L 

r 

, 2 S | 

^ 5 ) 

2 

. _ s 2 s 3 

f 

1 + 



ds = S 1 r- 

J 

0 

L 1 



L 3 L 2 


(2.42) 


(2.43) 


L 1 L 


= L + 


L 3L 


_r r L _ L 

— L — L h — — — 


(2.44) 



(2.45) 


Las funciones de forma referidas a un sistema coordenado local ubicado en el centra del 
elemento (Fig. 2.6), se obtiene reemplazando en las funciones de forma (2.22) y (2.23), 


x = x t + 




+q 


(2.46) 
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L2 



Fig. 2.6 Sistema coordenado local ubicado en el centra del elemento. 




X J -X 


L 


Xj-x, 


u. 


-q 


L 


L 

x j~ x i 2 q 

L L L 


1 q _ 1 q 

2 L~ 2 L 


(2.47) 


Nj (?) 


X-Xj 

L 


X : + 




v ^ j 


+ q~x t 


l - + «- 
2 L 


(2.48) 


Asl, las funciones de fonna, con respecto al sistema local, q , son: 


1 


"M-ri 


Que se pueden expresar en fiincion de x , 


VW = 
AM 


1 X 

2 L 

1 x 

- — l — 

2 L 


(2.49) 

(2.50) 


(2.51) 

(2.52) 


donde la variable dependiente varia de 
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EJEMPLO 2.3 

x j 

Evaluar la integral: j Nf(x)dx 

x i 

Para obtener la solucion se cambian los limites de integracion, lo que significa cambiar el 
origen de referencia de la funcion de forma /V. (a) — » N t (c/) , 

L L 

x i 2 2 

J N 2 (x) dx = J N 2 (<7) dq = J 



2 

2 

L 



2 

f 

i__<7 + <7l 

, 1 <? 2 
dq = — q + — ^ 

J 

L 

4 L L 2 

4 2Z 3Z 2 

2 




r l 


1 

"z 2 

Z 2 1 

1 

fz 3 

Z 3 1 


— 

— 

H 7 

1 

L2 

2 J 

2L 

4 

4 

3L 2 

8 

8 J 


Z Z 3 _L L _3L + L _4L _L 
4 12ZF~4 12~~12~~J2~J 


I_1 

2 Z 


dq 


(2.53) 


Resumiendo, las funciones de fomia del elemento linea son: 

«,(*) = 

W,(x) = . 


X - -x 

Z 

x-x. 


Z 


Sistema Local “s’ 




Sistema Local “q v 


lK ’ 2 L 

N-(x) = — + — 
]K ’ 2 L 


(2.54) 

(2.55) 

(2.56) 
(2-57) 

(2.58) 

(2.59) 
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Las integrates de los ejemplos 2.2 y 2.3, seran dificil de evaluar si: 


J N 2 (x)dx=\ 



|2 


dx 





Xj 


- 2 XjX + x 2 



L 2 


IXjX X 

xx 1 

1 2 3 


L 1 


: [x.-x ; ]-x.[x 2 -x, 2 ] + ^[; 



2.5 SISTEMAS COORDENADOS NATURALES 

Un sistema coordenado local en S y q sc puede convertir en un sistema coordenado 
natural. Un sistema coordenado natural, es un sistema local que permite la especificacion de 
un punto por un numero adimensional cuya magnitud absoluta nunca excede la unidad. 


L/2 L/2 



q 

1 



Fig. 2.7 Sistema coordenado local. 

Para el sistema con coordenada q que varia de a ^2 ^'8- 2.7), realizando el 

siguiente cambio de variable, 

f = y = ^ (2-60) 

2 

Se logra que los limites de integracion para la nueva coordenada £ sean de -1 a +1 (Fig. 

2.8). 
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C =-1 C =+1 


L/2 



c 

i 

r 


Fig. 2.8 Sistema coordenado natural. 


Por lo que la coordenada independiente es: 


q = 


4^ 


Reemplazando en la funcion de forma (2.49): 


N i (C) = --^ = --— = --^- = -{ 1-C) 

’ 2 L 2 2L 2 2 2 v ^ 


Reemplazando en la funcion de forma (2.50): 

& 

1 ? 1 2 1 Q 1 ^ 1 n ~ 

N.(C) = - + — = - + — = - + — = - + — = -(l + C) 

‘ 2 L 2 L 2 2Z 2 2 2 

Por lo tanto, las funciones de forma en funcion de coordenadas naturales son: 


(2.61) 

(2.62) 

(2.63) 

(2.64) 

(2.65) 



(2.66) 


(2.67) 


Para realizar estas integrales es necesario efectuar un cambio en los limites de integracion, 
de la siguiente manera: 


a 


2 Ci ? T +1 

J r{q)dq = \q{<;) 


(2.68) 


donde g(Q es r(q) escrita en terminos de £ . La ventaja de hacer el cambio de variables es 
que los limites de integracion son ahora de - 1 a + 1 . 
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Las matrices individuates de los elementos en ocasiones se integran en forma numerica 
usando algoritmos como el de Gauss-Legendre el cual integra en puntos con coeficientes 
definidos en el intervalo -1 a +1. 

Para un sistema de coordenadas S que varia de 0 a L, se pueden cambiar los limites por 
medio de relaciones de longitud (Fig. 2.9). 


11 12 





a 

s 

1 



j 


Fig. 2.9 Sistema coordenado local. 


Si S es la distancia del nodo i a un punto cualquiera del elemento, entonces las relaciones l\ 
y h estan definidas por: 


_ L - S _ a 
L ~ L 


(2.69) 


h= f (2-70) 

El par de coordenadas l\ y h estan relacionadas con la ecuacion: 

/j+/ 2 = 1 (2.71) 


Asi: 


h 


L-S S 
L ~ L 



que son identicas a las funciones de forma anteriormente definidas (2.38) y (2.39): 

W,(S) = 1-| = 7, 


(2.72) 

(2.73) 


(2.74) 

(2.75) 
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La utilidad de este sistema coordenado esta asociada con la evaluacion de integrales del 
tipo: 

L 

^N“{S)N b J {S)ds (2.76) 

0 

la cual involucra productos de funciones de forma. De la regia de cambio de variable y las 
relaciones: 



N,(S) = k 

(2.77) 


Nj(S) = 1 2 

(2.78) 


s=li 2 

(2.79) 

obtenida de (2.73): 

1 = s 

2 L 

(2.80) 

Diferenciando (2.79), 

dS = dl 2 L 

(2.81) 


_ dS 

J-j — 

dl 2 

(2.82) 

de (2.71): 

l\ + 1 2 = 1 

(2.83) 

Despejando /,, 

/j =1 -l 2 

(2.84) 

Sustituyendo en la integral (2.76), tenemos: 



L 1 

J (5) N) (5) ds = J l\l 2 Ldl 2 

0 0 

(2.85) 


1 

= Lj(l ~l 2 )“ l 2 dl 2 

0 

(2.86) 

Usando la formula de integracion: 



l r (z+w) 

(2.87) 

Donde, 

r(n + l) = n! 

(2.88) 
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Segun Abramowitz y Stegun (1964), tenemos: 


if rtdi -i r 0 +1 ) r ( 6+1 ) 

l 1 2 2 ~ r(a + b + l + l) 

(2.89) 

L f lyl 2 dl 2 = L a[b[ 

l («+* + !)! 

(2.90) 


Esta formula es importante para evaluar integrales que involucran longitudes de elementos 
y productos de sus potencias. 


EJEMPLO 2.4 
Evaluar la integral: 


Cambiando las variables de la integral, 


Como: 


entonces: 


| N f (x)r/x 

Xj 

(2.91) 

x i L 

\N?(x)dx=jN?(S)ds 

Xj 0 

(2.92) 

W,(S) = /, 

(2.93) 

n,(s)=i 2 

(2.94) 

L 1 

\N;(S)ds=L\l;lldl 2 

(2.95) 


0 0 


donde a = 2 y b = 0 . Empleando la formula de integracion (2.90), 


L J l[l\dl 2 
0 


_ r(2 + l)r(0 + l) _ 2!0! 

(2.96) 

r(2 + 0 + l + l) (2 + 0 + 1)! 

2-1 I , 

= L 

3-2 3 

(2.97) 


EJEMPLO 2.5 
Evaluar la integral: 

L 

^Nf(S)Nj(S)ds (2.98) 

0 
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Cambiando las variables de la integral, 


Jw/(S)w;(S)A 

0 


= lJ r;i';Ldi = l 

o 

_ L 3-2-2-1 
6-5-4-3-2 


3!2! 

(3 + 2 + 1)! 

L 

60 


(2.99) 

( 2 . 100 ) 


2.6 RESUMEN DE SISTEMAS COORDENADOS Y LIMITES DE INTEGRACION 
PARA UN ELEMENTO UNIDIMENSIONAL 


Tipo de sistema 

Coordenada 

Funcion de fonna 

Limites de 


variable 


integracion 

Global 

X 

ll 

J«S 

I> | 

X 

Xi ; xj 



Js; 

ii 

X 

1 

_x 


Local 

5 

N =\—— 

' L 

nJ- 

1 L 

0 ; L 

Local 

q 

N = — 

1 2 L 

N . = — + — 

J 2 L 

L L 

2 ’ + 2 





Natural 

<r 

+ 4 (1+f) 

-1 ; +1 



N t =h 


Natural 

h 

Nj=1 2 

0 ; 1 
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2.7 INTEGRACION NUMERICA EN UNA DIMENSION 

La evaluacion numerica de la integral: 

b 

I = J* / (x)<7x 

a 

es conocida como cuadratura. Las cuadraturas mas comunes en el MEF son 

a) Newton-Cotes 

b 

I = | f(x)dx 

a 

b) Gauss-Legendre 

i=\f(()d( 

-i 

2.7.1 CUADRATURA DE NEWTON-COTES 

Es un metodo numerico de integracion exacta para un polinomio de grado con una 
variable, donde n + 1 es el numero de funciones a evaluarse. El metodo se define como: 

/ = j/(x)<ix = £( Cn ,)/(*,) (2.104) 

a '=° 

Las abscisas x, son los puntos en los cuales la funcion es evaluada y se suponen igualmente 
espaciados (Fig. 2.10), Cn, son los pesos 6 funciones de peso. Supongamos que fix) es 
evaluada en 4 distintos puntos: Xo, X/, X 2 y Xj; espaciados a intervalos iguales, h . 

g (x) = Cj + C 2 x + C 3 x 2 + C 4 x 3 (2.105) 

La funcion g(x) se puede obtener por interpolacion Lagrangiana de fonna tal que, 

g( x ) = f(x 0 )L 0 (x) + f(x l )L l (x) + f(x 2 )L 2 (x) + f(x 3 )L 3 (x) (2.106) 

donde: 


o 

H 

i i 

~ * 

° i 

1 H 

A iu 

(2.107) 

L / A (X-X 0 )(X-X 2 )(X-X 3 ) 

(xi-x 0 )(xi-x 2 )(xi-x 3 ) 

(2.108) 


( 2 . 101 ) 


( 2 . 102 ) 


(2.103) 
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X„ X, x 2 x 3 

u h -1 


Fig. 2.10 Funcion de interpolacion. 


L / A (X-X 0 )(X-X 1 )(X-X 3 ) 

{x 2 -x 0 )(x 2 -x, )(x 2 -x 3 ) 

(2.109) 

L / A (x-X 0 )(x-X 1 )(x-X 2 ) 

( x 3 — X Q ) ( X 3 Xj ) ( x 3 — x 2 ) 

(2.110) 


La integration numerica puede realizarse usando interpolacion polinomica g(x) para cuatro 
puntos (Fig. 2.11). 



•N 

X X 


Fig. 2.11 Interpolacion polinomica. 


Considerando el espaciamiento entre puntos, 
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Sustituyendo en (2.104), 


donde: 


h 


( b-ci ) 

^3 


( 2 . 111 ) 


]f(x)dx = Y j \L, (. x)dx ( x ( ) + E 3 

a _ a 

= Z C 3 if( X i) + t 


- error 


( 2 . 112 ) 

(2.113) 



con i = 0, 1, 2, 3 


(2.114) 


= C 3o/(*0 ) + C 3 J (x, ) + C 32 /(x 2 ) + C 33 f (x 3 ) 


La formula general para las funciones de peso es: 


(2.115) 


C nl =hAW i (2.116) 

Los valores de A y de W t se obtienen de la Tabla 2.1. 


Tabla 2.1 Coeficientes de peso auxiliares para el calculo de cuadratura de Newton 


n 

A 

Wo 

W] 

w 2 

W 3 

W 4 

Error 

1 

Vi 

1 

1 




-1/2 

2 

1/3 

1 

4 

1 



-1/90 

3 

3/8 

1 

3 

3 

1 


-3/80 

4 

2/45 

7 

32 

12 

32 

7 

-8/945 


EJEMPLO 2.6 

Evaluar numericamente la integral: 

4 

/ = j" x~ dx (2.117) 

1 

cuya solucion exacta es 21 unidades cuadradas. 

a) Considerando n = 1 

b) Considerando n = 2 

c) Considerando n = 3 
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1. Considerando n = 1 (Fig. 2.12). 



Fig. 2.12 Interpolation para n = 1 . 


Calculo del espaciamiento, 


De la Tabla 2.1 se obtiene: 


, b — a 

h = 

n 


h 


4-1 

~r 


= 3 


A 


\_ 

2 


(2.118) 

(2.119) 


( 2 . 120 ) 


Pesos: 


Sustituyendo en (2.1 17), 


C o =hAW 0 = 3 


Cj = hA W x =3 




(1) = 1.5 
(1) = 1.5 


/ = C 0 /(1) + C i /(4) = 1.5(1) 2 +1.5(4) 2 


( 2 . 121 ) 

( 2 . 122 ) 

25.5 ^21 
(2.123) 


Esta es la regia trapezoidal la cual es valida solo para funciones lineales. 


2. Considerando n = 2 (Fig. 2.13). 

Espaciamiento entre puntos: 

h = ~~ = 1-5 (2.124) 

De la Tabla 2.1 se obtiene: 
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A = | (2.125) 



Fig. 2.13 Interpolacion para n = 2 . 


Funciones de peso: 


C 0 = hAW 0 = (l.5)f^-l(l) = 0.50 

V 

C,=hAW,={ l-5)ff|(4) = 2.0 
V-V 


1 \ 

C 2 = hAW 2 =( 1.5)| - (l) = 0.50 

3J 


Sustituyendo en (2.1 17), 


(2.126) 

(2.127) 

(2.128) 


/ = C 0 /(l) + C 1 /(2.5) + C 2 /(4) (2.129) 

/ = 0.5 (l) 2 +2(2.5) 2 +0.5(4) 2 =21 = 21 (2.130) 

Esta es la regia de Simpson 1/3, la cual es exacta para funciones cuadraticas y lineales. 

3. Considerando n = 3 (Fig. 2.14). 



Fig. 2.14 Interpolacion para n = 3 . 
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Espaciamiento entre puntos: 

h = ^Y~ = l-° (2.131) 

De la Tabla 2.1 se obtiene: 

A = - (2.132) 

8 

Las Funciones de peso son: W 0 = 1 , W l - 3 , W 2 - 3 , W 3 - 1 . Sustituyendo en: C n -hAW n . 

(2.133) 


Sustituyendo en (2.1 17), 


6 bien: 


^=w(D (, )=f 

c --w(i) (3 )-f 

c »=( , # ) =i 


I = C 0 /( 1) + CJ{ 2) + C 2 /( 3) + C 3 /(4) 


3 '-\2 , ^ , 9,-^2 , 3 ,.. 2 


(2.134) 

(2.135) 

(2.136) 

(2.137) 


I = -(1 y +-{2 y + o (3 y + o (4 y =21 = 21 

8 o o o 


(2.138) 

Esta es la regia de Simpson 3/8, la cual es valida para funciones cubicas o de orden inferior. 


2.7.2 CUADRATURA DE GAUSS -LEGENDRE 


Este es el procedimiento mas usado y popular en el MEF. Supongamos una funcion 
/ (^) para la cual se desea calcular la integral en el intervalo (-1, +1), es decir: 

+i 

/ = J/(^ (2.139) 

-i 

La regia de integracion o cuadratura de Gauss-Legendre expresa el valor de dicha integral 
como una suma de productos de los valores del integrando, en una serie de puntos 
conocidos en el interior del intervalo por unos coeficientes (pesos) determinados. Es decir, 
para una cuadratura de orden “p” se tiene: 
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I r = 'tf(C,W, (2.140) 

i=l 

Donde W; es el peso correspondiente al punto de integracion i y p es el numero de dichos 
puntos. La cuadratura de Gauss-Leg endre de orden p integra exactamente un polinomio de 
grado (2p - l) 6 menor. 

En la Tabla 2.2 se muestran las coord enadas £",• y los pesos W t para las seis primeras 
cuadraturas de Gauss-Leg endre. 

Observese que los puntos de integracion estan todos expresados en el espacio nonnalizado 
-l<£ <+l, el cual resulta de gran utilidad para el calculo de matrices de elementos 
referidos a coordenadas naturales. 

La popularidad de la cuadratura de Gauss-Legendre se debe a que utiliza un minimo 
numero de puntos de integracion para conseguir un error detenninado en el calculo, por 
consiguiente minimiza el numero de veces que hay que calcular el valor de la funcion a 
integrar. 

Tabla 2.2 Coeficientes de peso auxiliares para el calculo de cuadratura de Gauss 


p 

±Ci puntos de colocacion 


1 

0 

2 

2 

i+ 

1,1 



0 



9/8 


3 

±VO60 



5/9 , 5/9 



±0.33998 

10435 

84850 

0.65214 

51548 

62546 

4 

±0.86113 

63115 

94053 

0.34785 

48451 

37454 



0 


0.56888 

88888 

88889 

5 

±0.53846 

93101 

05683 

0.47862 

86704 

99366 


±0.90617 

98459 

38664 

0.23692 

68850 

56189 


±0.23861 

91860 

83197 

0.46791 

39345 

72691 

6 

±0.66120 

93864 

66265 

0.36076 

15730 

48139 


±0.93246 

95142 

03152 

0.17132 

44923 

79170 
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EJEMPLO 2.7 


+1 

Calcular la integral: / = j / (x)c(x , considerando p = 3 , y donde / (x) , es: 


-i 


m= 


X 

.2 


X 


+1 

I = \f{ x )dx=wj{c i )+iv 2 f{c 2 )+w 3 f{c 3 ) 


-1 


De la Tabla 2.2, para p = 3 , W x = ^,W 2 = ^ 


0 Sl/(x) = l 


ii) SI /(x) = x 


/ = -/(o)+-/(Vo^o)+-/(-Vo^) 


T 9 5 5 

/ = — i h — = 2.23 exacta 

8 9 9 


/ = -(VO60 ) + |(-V(X60 ) = 0 


exacta 


Hi) SI /(x) = x ? 


/ = ^-(0.60) + ^(0.60) = y exacta 


/v) Sl/(x) = . 


/ 


^(0.60)^ + ^(-0.60)^ =0 


exacta 


(2.141) 


(2.142) 


(2.143) 


(2.144) 


(2.145) 


(2.146) 
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v) Si f(x) = 


x 


vi) Si /(x) = x 5 


/ = |(0.60) 2 +|(-0.60) 2 =| exacta (2.147) 


vii) Si / (x) = x 6 


/ - ^ (0.60) 2 + ^ (-0.60) 2 =0 


exacta (2.148) 


/ = |(0.60) 3 +^(-0.60) 3 =0 


exacta (2.149) 


2.8 INTEGRACION NUMERICA EN DOS Y TRES DIMENSIONES 


En ocasiones es necesario calcular: 

i i 

I = \\f(C,rj)dCdri (2.150) 

- 1-1 

donde 2' y // pueden ser coordenadas cartesianas 6 curvilineas. La integral (2.150) se 
puede expresar, 

I = )^Hjf(Cj,fl)dV = )g{v)dV (2.151) 

-i j = i -i 

n m n 

I = £ H.g (n , ) = £ t Hfljf 7 , ) (2.152) 

(=1 i=\ j= 1 

Para el caso tridimensional: 

+ 1 + 1+1 

,S)d£,drj,dS (2.153) 

-l-i-i 

que se puede expresar, 

l m n 

I = I II H.HjHJ ((„ nj 4,) (2.154) 

i=l 7=1 k = 1 

En muchos casos no hay ventaja en usar diferentes puntos de integracion, si no que se 
supone que / = m = n . Si la funcion es integrada con / = m = n = 3 la suma equivale a 
3x3x3 puntos, la cual es exacta para polinomios de quinto orden en cada direccion. 
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Altemativamente, la integral se puede escribir: 

+ 1 +1 +1 Yl 

I = \\\f(C,^)d^di 1 ,d^ = ^W i f{C i ,i li ^ i ) (2.155) 

-i -i -i <=i 


EJEMPLO 2.8 

Considerando n = 2, desarrollar la formula de integracion y evaluarla para el elemento 
cuadrilatero de la figura 2.15. Los valores de los pesos se obtienen de la Tabla 2.2. 


+1+1 


^2x2 \ \ f{C,r?)dC dr/ 

- 1-1 

+iT+i 

| \f((,r,)d( d n 


/ = 


1 = \ [ W J {Cl, v) + w 2 f {Cl, r/)yr/ 


(2.156) 

(2.157) 

(2.158) 


Los pesos son constantes y pueden salir del operador integral, 

I = W x \WJ^ x ,r 1 ^W 1 f{^r 1 ^\ 

+W 2 \WJ^ 2 ,7 lx ) + W 2 f^ 2 ,7 l2 )\ 


(2.159) 

Asi, la integral queda definida por, 

^2=^ 2 /fe.7,) + ^/(^7 2 ) + ^2/fc.7l) + ^2 2 /(^.7 2 ) 

(2.160) 

de la Tabla 2.2, tomando los pesos y los puntos de evaluacion, para n=2: 

W l = 1 771 = “73 (2 ' 161a) 

w 2 = i Vl=+ ifs (2 ' 161b) 

Es decir, la integral sera evaluada en el elemento cuadrilatero de la figura 2.14, en los 
puntos 4), >h y ^ 2 , r h 
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’ h “ + V3 










V3 


Fig. 2.15 Elemento cuadrilatero de integracion. 


Sustituyendo los valores de los pesos y los puntos de evaluacion en: 

A rf =^ 2 /fe.7i) + ^/(^%) + ^/fe.7i) + ^2 2 /fc.7 2 ) (2-162) 

Tenemos: 


1 2,2 = CD 


f 1 \ 


1 


V -V 


f i \ 


+ 


1 


V 


+(1)(1) 


f 1 \ 


1 


v 


+ 


6 1 ^ 

H 

v 3 j 


+ 




( i A 

+ - 


( i > 


+ 


v V3 J l V3 


v VJ y 


+ 


( 1 ) 


A i '\ 

+ - 


+ 


( i 
+ - 


v S) l V3 


v VJ y 


Realizando operaciones, 


I = - 0.244016 + 0.910684 - 0.244016 + 0.910684 


(2.163) 

(2.164) 


el valor exacto es 


7 = 1.333 



(2.165) 

(2.166) 
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EJEMPLO 2.9 


Considerando n = 3 desarrollar la formula de integration: 




3x3 




- 1-1 


F7 


(2.167) 


+ + (2.168) 

-1 

= w l [wj{c l ,n l ) + wj{^ + wj{c l ,? 7i )] 
+w^wj(c 2 ,n i ) + wj{C2^2)+wj(^ 2 ,n i )] 

+W 3 [Wj{^ 3 ,r ll ) + W 2 f(C^2) + W 3 f{C 3 ,rj 3 )] (2.169) 
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CAPITULO III 

ELEMENTOS UNIDIMEN SION ALE S Y BIDIMEN SION ALES 

3.1 ELEMENTO FINITO BARRA 

Este elemento tiene solo deformaciones axiales. 


^11 0 


AE 

dx 



Fig. 3.1 Deformaciones axiales. 


Considerando que A es el area de cortante, E el modulo de elasticidad del material y u j los 
desplazamientos del nodo “i” en la direccion x. Tomando un elemento diferencial: 


F = Act, 


F = Act . + A 


ax 


Fig. 3.2 Elemento diferencial. 

donde o x es el esfuerzo normal en la direccion x. Del equilibrio, 

2>,=o 

-Act .+ Act. + A = 0 


dx 


A^ = 0 
dx 

De la ley de Hooke 

<*x = E 

donde la deformacion unitaria en la direccion x, es: 


e 


du 

dx 


d cr x 
dx 


(3.1) 

(3.2) 

(3.3) 


Sustituyendo en (3.2), 

^ du 
cr =E — 
dx 

y sustituyendo en (3.1), 


(3.4) 
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AE 


d 2 u 

dx 2 


= 0 


(3.5) 


Que es la ecuacion diferencial de una barra sometida a esfuerzo axial. 

La ecuacion diferencial se resolvera por el Metodo de Residuos Pesados, particulannente 
por medio del metodo de Galerkin. 


1 . Se supone una aproximacion a la solucion: 


u* = + N 2 Q> 2 (3.6) 

donde <f>j y <P 2 son los valores nodales. La ecuacion (3.6) tambien se pude escribir, 

n 

u = N\U X + N 2 u 2 = ^ N i u i (3.7) 

i=i 

donde las funciones de interpolacion o de forma, referidas a un sistema coordenado local, 
son: 


7V,= 


!-* 


N =1 

2 L 


Sustituyendo en la aproximacion de la solucion (3.7), 


u = 


"l-* 

W| 

+ 

X 

L_ 



_L_ 


2. Sustituimos en la ecuacion diferencial (3.5), 

, d 2 u 


AE- 


dx" 


= R 


(3.8) 

(3.9) 

(310) 

(3.11) 


3. Se pesa el error ^RW^Q. = 0 

Q 


4. Por el metodo de Galerkin: 



* 

d u 
dx 2 


\W : dx = 0 


los pesos son las funciones de fonna, que sustituyendo en (3.12), 

' >2 * 
d u 


EAj 


dx~ 


Nclx = 0 


De la integracion por partes: 


d[A-B] = AdB + BdA 
Despejando: BdA = d[A- B]- AdB 
Integrando : J BdA = A- B - 1 AdB 

dN. 

con: B = N t dx dB = — '-dx 

dx 


dA = 


d 2 u 

dx 2 


A 


du 

dx 


(3.12) 

(3.13) 

(3.14) 


(3.15) 
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Empleando la formula (3.15), para resolver (3.14), 

du 


L j2 


EA f N dx = EA 

J dx 2 1 


N t 
dx 


L r du dN. 


0 ' dx dx 


L dx 


en donde: 


a = Ndx 


A dN ‘ A 

du = dx 

dx 


, d 2 u , du 

dv = — —dx v — 

dx dx 


De (3.16), 


Despejando la integral definida, 


EA 


du 

dx 


N, 


I ^ 

0 J 


du dNf 
dx dx 


dx 


= 0 


E A\E d JL dx=EA d JU N . 


Como: 

Sustituyendo en (3.18), 


u = N ] u ] + N 2 u 2 ; 


dx dx 

y 


dx 

du* 
dx 

dN : 


= 0 


EA f — [N j u i + N 2 u 2 ] — L dx = EA [e r ] N t 
{dx “ dx 


como 


cr -Eg, tenemos: 


como: F = Aa x 

matricialmente 
como: i=l,2 

Para i = 1, en (3.21), 

Para / = 2, en (3.21), 


EA\- 


f dN, 

dN. dN, 

dx = Act N, 

-U l H -U 2 

- dx 

dx dx 

X l 


ea\ 


dN : 


dx 


dN x dN 2 


dx dx 


dx = FN: 


EA 


f dN\ 

~dN x 

dN 2 " 

Uj 

<: 

ii 

"is 



• dx 

dx 

dx 

u 2 



r dN, 

~ dN x 

dN, 1 

u \ 

£ 

II 

1 2 

2 


l dx 

dx 

dx 

_ U 2 _ 

z. | 


Las ecuaciones (3.22) y (3.23) matricialmente se pueden escribir, 


(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 

(3.22) 

(3.23) 
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~ dN x “ 




L 

dx 

dN l dN 1 

u \ 


~n x 1 

EA\ 




fen 

II 

3 


J 

0 

dN 2 

dx dx 

u 2 


n 2 \ 


dx 





Reorganizando: 



~ dN x ' 







L 

f 

dx 

[ea] 

^3 

1 

1 

CN 

£ 

^3 

u \ 

& 

ll 

>r 

L 

J 

0 

dN 2 


dx 

dx 

u 2 


n 2 

0 


dx 








Que tiene la forma: 

donde la matriz elemental de rigidez es: 

De (3.26) se puede ver que, 

como las funciones de forma (3.8), (3.9): 




1 -* 

L 


N =- 
2 L 


sus derivadas, 

sustituyendo en (3.25), se tiene: 

EA\ 

realizando operaciones: 

ea\ 

0 

como L (longitud del elemento) es constante: 

Ea 
L 2 


dN 2 _ 1 
dx L dx L 



1 “ 




L 

r 

L 

1 

1 

1 \u, ] _ ( 

f 


+ — 

K 

ll 

§ 

j 

0 

1 

+ — 

L 

L 

1 K J l 


L 





Ni(L) ~N r (0) 
7V 2 (Z) -N 2 { 0) 


f i 

1 ■) 

! 


L 2 

L 2 

II 

, ^ 

( 0 -f 

1 

1 

[« 2 J 1 

vl 0 y 


Tf , 




j( 1 -1) 

f ll, 1 

'-f 

f 

5 

II 


J ol-l 0 

KJ 

+i 


(3.24) 


(3.25) 


J[ B f[ D ][ S ]{ u pv = {P} 

V 

(3.26) 

[r] = J[5f[D][5]iv 

V 

(3.27) 

MH=M 

(3.28) 


(3.29) 


(3.30) 


(3.31) 


Como los desplazamientos son constantes en toda la longitud de la barra, y llamando: 


F x =-F 

F 2 =F 
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Ft 


Fj 


x, u 


Fig. 3.3 Fuerzas en la barra. 


La expresion (3.31) se puede escribir, 

EA 

IF 

6 bien, en coordenadas locales: 


F -l) 

\ u \ 1 r 

r^i 


( \\dx- 


v-1 K 

Mo 

A. 


EA 

( 1 

-1) 

h 1 = 

>r 

L 


K 

KJ 

A. 


Donde la matriz elemental de rigideces es, 




EA 

L 


1 -1 

-1 1 


(3.32) 


(3.33) 


(3.34) 


[Referencia] ‘Analisis Avanzado de Estructuras’, apuntes del Ing. Hugo Hernandez Barrios. 


EJEMPLO 3.1 

Determinar los desplazamientos en los nodos, las defonnaciones y los esfuerzos para los 
elementos del miembro mostrado en la figura. Considere el modulo de elasticidad, E, y el 
area, A, de la seccion transversal de los elementos, como constantes. 



Fig. 3.4 Elementos y nodos de una viga en cantiliver. 


Anteriormente se determino que la matriz de rigidez elemental, para el elemento 1 , 

u ] u 2 


para el elemento2, 


EA 

L 




u \ 
ll 2 



EA 

L 


-1 


« 3 

-n 


U-, 

u 3 


La matriz de rigidez total: 
es decir: 
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De la relation: 


u, ll-, u. 


' 1 -1 0 '\u l 

— 1 1 + 1 — 1 ^2 
0 -1 1 u. 




F, (l -1 0 1 m, 

r-, EA 

F 1 = -1 2 -1 < u 1 

L 2 

F \ 0 -1 1 u , 


Las condiciones de frontera son: 

Mj =0 en el apoyo que esta empotrado. 
F, = P en el cantiliver. 


1 -1 010 


0 = — -1 2 -1 \\u. 


0 -1 1 \u 


Resolviendo el sistema de ecuaciones anteriores: 

„ EAr , 

P = —[-u 2 +u 3 \ 


de donde: 


0 = 2 u 2 - u 3 
u 3 = 2m ? 


„ EA , ^ i EA 

P = [-m 2 + 2r< 2 J = -j-u 2 ; 


despejando los desplazamientos nodales: 


i PL 
u , = 2 


. . r, E A EA PL 

para la nierza : P = m, = = -P , 

1 L ~ L EA 


F\ = ~P 


DEF ORMACIONES ELEMENTALES: 


Elemento 1: 




1 ll «i 


L L u 


L L 
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ESFUERZOS ELEMENT ALES: 
Elemento 1: 

Elemento 2: 




3.1.1 MATRIZ DE RIGIDEZ EN COORDENADAS GLOB ALES 


La matriz de rigidez de un elemento barra en coordenadas locales es: 

M=-f‘ ‘‘1 

L J L -1 1 


(3.35) 


Cuando el elemento esta inclinado un angulo a , es necesario calcular la matriz de rigidez 
en coordenadas globales. 



Fig. 3.5 Relacion entre desplazamientos. 
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Sean x -y los ejes globales y x-y los ejes locales de la barra. Se puede ver que la 
relacion entre los desplazamientos u y v y los desplazamientos globales u y v es: 

u = u cos a + vsena 
v = v cos a-u sena 


que se pueden escribir: 

r cos a sena 
-sena cos a 

Asl, para los dos puntos nodales, dicha relacion se puede escribir: 


(3.36) 


6 bien: 


u l 


cosa 

sena 

0 

0 


u x 

V 1 


-sena 

cos a 

0 

0 


v 'i 


> = 





< 


u 2 


0 

0 

cos a 

sena 


u 2 

V 2. 


0 

0 

-sena 

cos a 


V 2, 


[u] = [T]{u e ] (3.38) 


En el caso de que unicamente sea necesario calcular los desplazamientos en la direccion 
del eje de la barra (desplazamientos u ] y u 2 (ejes locales), la ecuacion (3.37) se puede 
escribir, 



\ u x 

L_ 

cosar 

sena 

0 

0 

i U 2. 


0 

0 

cos a 

sena 


V 2. 


(3.39) 


donde u [ , vj , u 2 y vj son los grados de liberad de la barra en coordenadas globales. 


Para el calculo de las fuerzas nodales en coordenadas locales, segun la siguiente figura, se 
tienen: 



Fig. 3.6 Fuerzas nodales en coordenadas locales. 
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>*]: 


X 


‘ 1 

0 

-1 

0 


_ EA 

0 

0 

0 

0 

f 2x 

~ L 

-1 

0 

1 

0 

F 2y 


0 

0 

0 

0 


(3.40) 


donde F lx , F ly son las fuerzas en direccion x y y, en el nodo 1, F 2x , F ly son las fuerzas en 
la direccion x y y, en el nodo 2. 

Si unicamente interesa calcular las fuerzas en la barra en la direccion del eje local, x: 

(3.41) 

Por tanto, las fuerzas nodales en la direccion local, x, de la barra, se obtienen de sustituir la 
ecuacion (3.39) en la (3.41), 

r ■ 
u. 


(3.42) 


X" 

EAl 

r i -n 

j“d 

x_ 

- L 

x 0 

KJ 


X" 

_EA\ 

( 1 

-n 

cos a 

sena 

0 

0 

1 

K) J 

5-! 

1 

~ L 

v-1 

K 

0 

0 

cos a 

sena 


donde u\ , vj , u 2 y v' 2 son los grados de liberad de la barra en coordenadas globales. 
Se puede demostrar que la matriz de rigidez en coordenadas globales es: 

M=[ r fM[ 3 ’] 


(3.43) 



c 

cs -c 

-cs 

EA 

cs 

s 2 -cs 

-s 2 

L 

- c 2 

-cs c 2 

cs 


-cs 

-s 2 cs 

s 2 

: L-- 

= J(x 2 

-X { ] + [y' 2 

-y\. 


c = 


s = 


L 

(y'i-y\) 

L 


= cos or 


= sena 


(3.44) 


EJEMPLO 3.2 

Calcular los elementos mecanicos en las barras de la annadura plana de la figura, cuyos 
elementos tienen las siguientes propiedades: 


Elemento 

Area (in 2 ) 

Longitud (in) 

(kips/in) 

cosa 

sena 

1 

12 

240 

1500.0 

0.800 

0.600 

2 

12 

180 

2000.0 

0.600 

-0.800 

3 

8 

300 

800.0 

1.000 

0.000 
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Fig. 3.8 Topologia del modelo. 


Empleando la matriz (3.44) formamos la matriz de rigidez elemental referida a ejes 
globales, para cada barra, 


BARRA 1 


BARRA 2 


1 

2 

5 

6 



M, 


u 3 

v 3 



960 

720 

-960 

-720" 

U \ 


720 

540 

-720 

-540 

v \ 


-960 

-720 

960 

720 

u 3 


-720 

-540 

720 

540 

V 3 


5 

6 

3 

4 



u 2 

V 3 

u 2 

V 2 



720 

-960 

-720 

960 " 

u' 

-960 

1280 

960 

-1280 

v 3 

-720 

960 

720 

-960 

u\ 

960 

-1280 

-960 

1280 

v 2 


kips/in 


kips/in 
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BARRA 3 


4 



«i v i 

800 0 

0 0 

-800 0 

0 0 


U 2 V 2 


-800 0 

0 0 

800 0 

0 0 


Mj 


La matriz total, referida a coordenadas globales es: 


1 

2 

3 

4 

5 

6 


M, 

v i 

u 2 

v 2 

u 3 

V3 


1760 

720 

-800 

0 

-960 

— 720 

Mj 

720 

540 

0 

0 

-720 

— 540 

v| 

-800 

0 

1520 

-960 

-720 

960 

U 1 

0 

0 

-960 

1280 

960 

-1280 

*2 

-960 

-720 

-720 

960 

1680 

-240 

u 3 

-720 

-540 

960 

-1280 

-240 

1820 

V 3 


De la relation: 

Mapp'd 

Por las condiciones de frontera o de restriction: 

En el apoyo fijo: 

m, = 0 

v, = 0 


En el apoyo movil: 


= 0 


Tenemos que la operacion: 


1760 

720 

-800 

720 

540 

0 

-800 

0 

1520 

0 

0 

-960 

-960 

-720 

-720 

-720 

-540 

960 


0 

-960 

-720" 


0 


R lx ~ 

0 

-720 

-540 


0 


R ly 

-960 

-720 

960 


u 2 


0 

1280 

960 

-1280 


0 

=■ 

R 2 > 

960 

1680 

-240 


u 3 


0 

-1280 

-240 

1820 


V 3. 


-30 


Condensando la matriz, 


1520 

-720 

960 


u 2 


0 

-720 

1680 

-240 

< 

u 3 

> = « 

0 

960 

-240 

1820 


vi 


-30 


kips/in 


kips/in 


Resolviendo para los grados de libertad en coordenadas globales, es decir para los 
desplazamientos nodales, 
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[U G ] = 


0.00000 

0.00000 

0.01800 

0.00000 

0.00408 

-0.02544 


in 


Sustituyendo en la relacion previamente deducida, 


"1760 

720 

-800 

0 

-960 

-720 ' 


' 0.00000 


X" 

720 

540 

0 

0 

-720 

-540 


0.00000 


R ly 

-800 

0 

1520 

-960 

-720 

960 


0.01800 


0 

0 

0 

-960 

1280 

960 

-1280 

5 

0.00000 

y — 

*2, 

-960 

-720 

-720 

960 

1680 

-240 


0.00408 


0 

-720 

-540 

960 

-1280 

-240 

1820 


-0.02544 


-30 


Resolviendo para el vector de cargas nodales, 


K 


0.000' 

R U 


10.80 

0 

> = < 

0.000 

R 7 

2 y 


19.20 

0 


0.000 

-30 

v ^ 


-30.0 

C J 


kips 


Para encontrar los elementos mecanicos en las barras, es decir, las fuerzas axiales (con 
referencia a ejes locales), se sustituye en (3.42), los valores de los desplazamientos nodales 
en cada barra, 


BARRA 1 


(3.42) 


X' 

EA\ 

f 1 

-n 

cos a 

sena 

0 

0 

i 

X 

1 

~ L 

v-1 

K 

0 

0 

COS1Z 

sena 


X" 

= 1500 

f 1 

-n 

"0.8 

0.6 

0 

0 ' 

x_ 


v-1 


0 

0 

0.8 

0.6 


0.00000 

0.00000 

0.00408 

-0.02544 


X" 


" 18.00 ' 



-18.00 


kips 


Con respecto a los ejes locales la Barra 1 esta en compresion. 


67 



FUNDAMENTOS BASICOS DEL ELEMENTO F1N1TO Hernandez Barrios Hugo FESA-UNAM 


BARRA 2 





' 0.00408 ' 

1X1 


r i -r 

"0.6 -0.8 0 0 


-0.02544 

j X 

= 2000 



< 

> 

x_ 


l-i K 

0 0 0.6 -0.8 


0.018000 




0.000000 


X" 


■ 24.00 “ 

x_ 


-24.00 


u 3 

V 3 

u 2 

V 2 


Con respecto a los ejes locales la Barra 2 esta en compresion. 


BARRA 3 


X" 

= 800 

f 1 

-n 

"1.0 

0 

0 

o' 

X. 


v-1 

h 

0 

0 

1.0 

0 


0.00000 

0.00000 

0.018000 

0.000000 




“-14.40“ 

X_ 


14.40 


Con respecto a los ejes locales la Barra 3 esta en tension. 


kips 


kips 


A continuacion se explica paso a paso como realizar el ejemplo de la annadura en el 
programa SAP 2000 V.10. 

ARMADURA DEL EJEMPLO 3.2 [P.65.] 

1.- CREACION DE UN ARCHIVO NUEVO 

File > New Model 
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Initialize Model from Defaults with Units [seleccionar unidades] 
Kip, in, F [kips, inches, farenheit] 
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Dar clic en 2D Trusses [armaduras] 



Number of Divisions 1 [numero de divisiones] 
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Height 144 [in] [altura] 

Division length 300 [in] [longitud] 



Palomear la casilla: Use Custom Grid Spacing and Locate Origin y dar clic en Edit Grid 


'Q\sM 


II » fw i I « 

X 


'a- r*. I PTTM I 


i%; Define Grid Data 




System Name 

|CSYS1 



Grid ID 

Ordinate 

Line T ype 

Visibility 

Bubble Loc. 

1 

xl 

•150. 

Primary 

Show 

End 

2 

x2 

0. 

Primary 

Show 

End 

3 

x3 

150. 

Primary 

Show 

End 

4 






5 






6 






7 






8 








Grid ID 

Ordinate 

Line Type 

Visibility 

Bubble Loc. 

1 

yl 

0. 

Primary 

Show 

Start 

2 






3 






4 






5 






6 






7 






8 








Grid ID 

Ordinate 

Line Type 

Visibility 

Bubble Loc. 

• 

1 

zl 

0. 

Primary 

Show 

End 

j 

2 

z2 

144. 

Primary 

Show 

End 


3 







4 







5 







6 







7 







8 








| Kip, in. F ~^1 


Bubble Size |60. 


JnrrM’ U 



Display Grids as 
(* Ordinates C Spacing 


r Hide All Grid Lines 
V Glue to Grid Lines 


Reset to Default Color 


Reorder Ordinates 


Locate System Origin... 


BE: 


Cancel | 
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En la segunda columna “Ordinate” cambiar xl a -192 y x3 a 108 


□00 


□ at 


fg. I 


•f*;. Define Grid Data 


si I n rrW 


|CSYS1 


| Kip, in. F ~^1 


X Grid Data 








Grid ID 

Ordinate 

Line Type 

Visibility 

Bubble Loc. 

Grid Color 

1 

xl 

-192 

Primary 

Show 

End 

2 

x2 

0. 

Primary 

Show 

End 

Wr 

3 

x3 

108 

Primary 

Show 

End 


4 







5 







6 







7 







8 







Y Grid Data 








Grid ID 

Ordinate 

Line Type 

Visibility 

Bubble Loc. 

Grid Color 

1 

yl 

0. 

Primary 

Show 

Start 

2 







3 







4 







5 







6 







7 







8 







Z Grid Data 








Grid ID 

Ordinate 

Line Type 

Visibility 

Bubble Loc. 

- 

1 

zl 

0. 

Primary 

Show 

End 


2 

z2 

144. 

Primary 

Show 

End 


3 







4 







5 







6 







7 







8 






d 



Display Grids as 
(* Ordinates C Spacing 


|“ Hide All Grid Lines 
f" Glue to Grid Lines 


Bubble Size 1 60. 


Reset to Default Color 


Reorder Ordinates 


Locate System Origin... 


Cancel | 


3510 


Palomear las casillas: Hide All Grid Lines y Glue to Grid Lines [derecha] 


M0 


QtS 

llB^. t l« 


M 


at. ft. ja-ia-iJMJ. 


•«. Define Grid Data 


0 


hhH’U 


|CSYS1 



Grid ID 

Ordinate 

Line Type 

Visibility 

Bubble Loc. 

1 

xl 

-192 

Primary 

Show 

End 

2 

x2 

0. 

Primary 

Show 

End 

3 

x3 

108 

Primary 

Show 

End 

4 






5 






6 






7 






8 








Grid ID 

Ordinate 

Line Type 

Visibility 

Bubble Loc. 

1 

yi 

0. 

Primary 

Show 

Start 

2 






3 






4 






5 






6 






7 






8 








Grid ID 

Ordinate 

Line Type 

Visibility 

Bubble Loc. 

- 

1 

zl 

0. 

Primary 

Show 

End 

j 

2 

z2 

144. 

Primary 

Show 

End 


3 







4 







5 







6 







7 







8 






d 


| Kip, in. F ~^1 



Display Grids as 
(+■ Ordinates C Spacing 

[7 Hide All Grid Lines 
[7 iGlue to Grid Lines! 


Reset to Default Color 


Reorder Ordinates 


Locate System Origin... 


Cancel | 


][n)B 
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Ok > Ok 



Escoger un lado para trabajar, en este ejemplo trabajaremos del lado derecho [dar clic en 
ventana derecha] 

2.- DEFINIR MATERIALES Define > Materials 


'M SAP2000 vl 0.0.1 Advanced -(Untitled) 

File Edit View Define Bridge Draw Select Assign 

H H fci 




Analyze Disglay Design Options Help 

*_ m . 
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Steel > Modify/Show Material 



Modulus of Elasticity >30 000 > Ok > Ok 
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3.- DEFINIR SECCIONES Define > Frame Sections 




Choose Property Type to Add > Add Circle 
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Choose Property Type to Add 


Find this property: 


l/Wide 


Add Circle 


W8X10 

W8K24 

W8X31 

WECK40 

W8X58 

W8X67 

W1CCX12 

W10X49 

W10X68 

W1CK88 

W10X112 

W12X14 


Add Box/Tube 
Add Pipe 
Add Rectan gular 


Add BU I Cover Plate 
Add PC Cone I Girder 
Add General 
Add Cold Formed C 


[global 


AP200C vIO 

0.1 Advanced - (Unfitted) 


□00 

File Edit View 

Define Bridge Draw Select Assign Analyze Dismay Design Options Help 



□ u H 

■n } ® & & & £> a 3-d »y xz yz < » *(• O -0 ?h 0 

'a . 

I • a . nh«- . 

life* l« 

V - " 1 i- s - 




X-Z Plane @Y=0 


300 


Z Plane ® Y=0 


Frame Properties 


Click to > Add New Property 



Section Name > Elemento 1 y 2 [nombre de la seccion] 
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Material > Steel [material] 
Diameter > 3.91 [diametro] 



Ok 



Properties 


Choose Property Type to Add 


| Import I /Wide Flange 


Find this property: 


|ELEMENT01Y2 


ELEMENT01Y2 


W8X10 

WEK24 

WEK31 

WSK40 

WEX58 

W8XG7 

W10X12 

W1CK49 

W1CK68 

W1QX88 

W1 0X112 


Id New Proper 


Add Copy of Property... | 


M odify/S how Property. . . | 


M APJOOOvlO.O • Advanced Untitled; 

File Edit View Define Bridge Draw Select Assign Analyze Dismay Design Options Help 




□ b? H ® | } S > P P & & & B3 1 3<I « iz . a- M I -o- - 0 - [ gj a | y „ . 

l | a ,i.| || £ ei h m y | -- ? | i, & ii ,| 


3-D Viev. 


Kx Z Plane* 


□00. 


Frame Properties 


Delete Property 


B GLOBAL j]|Kipin.F jj 


XO.OO YO.OO ZO.OO 
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Choose Property Type to Add > Add Circle 


% 

SAP2000 vIC 

0.1 Advanced - (Untitled) 

□E0 

File 

Edit View 

Define Bridge Draw Select Assign Analyze Dismay Design Options Help 

1 

□ & a 

«o $ (x) ► P> $£) £) $5 HI 3 ’d xy xz yz > Ct O O ?h 0 X - 

I -H , nnM' . 


« '' ■ ■ ~ 1 - & - - 

& 

ffl D V 

QEB : S5 x - z plans 1 v "0 

□00 


< 


X-Z Plane @Y=0 



X0.00 Y0.00 Z0.00 


Section Name > Elemento 3 [nombre de la seccion] 
Material > Steel [material] 


Global j t „ F d 


% 

SAP2000v1 0.0.1 Advanced (Untitled) 

□B0 

File 

Edit View Define Bridge Draw Select Assign Analyze Dismay Design Options Help 


Q BS H ® 1 « ; S ’ H 3-d »y yz dtf -0 O H * . 

I • a . nhM- . 

1 

. IjUu-M id ^ - 


w 

X D View -1001 SS ■ Z Plane ® Y-0 




X-Z Plane @Y=0 


X0.00 Y0.00 Z0.00 


1.L0BAL 1 
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Ok > Ok 



4.- ASIGNAR CONDICIONES DE APOYO, TIPOS DE SECCIONES Y CARGAS 


4.1 Asignar Condiciones de Apoyo 
Seleccionar el primer apoyo con el mouse 

H< SAP2000v10.0.1 Advanced -(Untitled) 

File Edit View gefine Bridge Draw Select Assign Analyze Dismay Design Options Help 

D G» at* “O } & ' P @ @ @ & H M * xz yz ... Cf to I & O I B | ,| jA 

M iu * SS »I . H it m Lti ¥ | •*« V, | i_ SE (2 . 
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Assign > Joint > Restraints 




Apoyo fijo > Dar clic > Ok 
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Seleccionar el segundo apoyo con el mouse 


* 

SAP2000v1 0.0.1 

Advanced - (Untitled) 


- US 

File 

Edit View Define Bridge Draw Select Assign Analyze Display Design Options Help 

— 

a a h Sf « 

(T ® ► P P P P P 1 

- & at | « o ?n 0 >6 . i-a- 

. nrrH - . 

! 

t « -i . 




is 

M D 

□00 

Joint Restraints 

“ ESS 



Assign > Joint > Restraints 



Apoyo movil > Dar clic > Ok 
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4.2.- Asignar Tipos de Secciones 
Seleccionar elemento 1 y 2 con el mouse 



Assign > Frame/Cable/Tendon > Frame Sections 


M SAP2000 vl 0.0.1 Advanced -(Untitled) 

File Edit View Define Bridge Draw Select Assign Analyze Display Design Options Help 



| GLOBAL £]| Kip, in.F 3 


X-194.12 Y0.00 Z283.G1 
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Seleccionar Elemento 1 y 2 > Ok 


'M SAP2000 vl 0.0.1 Advanced -(Untitled) 




File Edit View Define Bridge Draw Select Assign Analyze Display Design Options Help 


DBS 9 IS 


► P P P P P H 3-d xy Hz yi ng Cf \ 3} O | 0 | 7 


- 1 ,| || n fr H - 1 .| 


5 

-% 


\ 

\ 

* 

<7 

□ 

0 

,i» 

PS* 

clr* 

-> 


K 3 D View 





1%, Frame Sections 


= m 



X-200.52 Y0.00 Z57.30 |gL 0BAL ^]|Kip. in. F _^J 
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Seleccionar elemento 3 con el mouse 



Assign > Frame/Cable/Tendon > Frame Sections 


M SAP2000 vIO.O • Advanced (Untitled- 

BBS: 

File Edit View Define Bridge Draw Select Assign Analyze Dismay Design Options Help 

1 


DtJ HSSj|«; : S 

I Si . ... ' I : , ii 8H ► 


is) & & & H] 3-d xy xz yz Cf (ftf | O | 0 | . 

I -£* X, ±. g£! M= - 


nhH’U 


M 3-D •' 


:: X 


0B. 



Global 


X0.00 Y0.00 Z0.00 
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Seleccionar Elemento 3 > Ok 



4.3.- Asignar Cargas 

Seleccionar la articulacion donde sera aplicada la carga 


* 

d 

o 

> 

o 

© 

© 

< 

Advanced - (Untitled) 


BBB 

File 

Edit View Define Bridge Draw Select Assign Analyze Disglay Design Options Help 


□ u h a* « 

} & ► P P P P P P 05 3-d »y a yz 

o a? o o ?h 0 y a , 

|| I ’ □ . nrrW’ , 

Jl* 

5 a «i . 





M D 

CHS 

M* Frame Sections 

. sa 
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Assign > Joint Loads > Forces 



1 Points Selected 


X0.00 Y0.00 Z0.00 | GLOBAL | Kip, in. F 3 



Force Global Z > -30 > Ok [fuerza en el eje Z global] 
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5.- ANALIZAR [EJECUTAR] EL PROGRAMA 

Analize > Run Analisis 
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M AP2000 v1 0.0.1 A 

Jvanced - (Untitled) 


□ 00 : 

File Edit View Define 

Bridge Draw Select Assign Analyze Dismay Design Options Help 


1 

d & a & « 

} fi > J9 & $ & & H 3-d nil a yz Ct U « -0 ?,J 0 

'o . 

I • a . nh«- . 

- 

\L - M V • ' i. - 




* 



X-Z Plane @Y=0 


X0.00 Y0.00 Z0.00 


J GLOBAL ~3F 3 


Run Now > Click 

Save Model File As [crear un folder, poner nombre y guardarlo ahi] 
Save 



X-Z Plane @ Y=0 X0.00 Y0.00 Z0.00 “3F J 
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id 3 

Goode Trooa Relaaor 


a 

Antivrus_ES SAP 2000 8 
Nonlinear 

II 0 


Adobe Reader RecordNow! 


Norton My Network 
SystemWorks Places 


M SAP Analysis Monitor 

File N ame: C: \D ocuments and S ettings\R olando G arciaSMy D ocuments\S ap\Armadura 3. 2ST ulorial 3. 2\tutorial3. 2. sdb 
Start Time: 5/3/2007 1:48:15AM Elapsed Time: 00:00:03 

Finish Time: 5/3/2007 1 :48:1 8 AM Run Status: Done - Analysis Complete 






A 

LINEAR EQUATION SOL 

U T I O N 


01:48:15 


FORMING STIFFNESS AT ZERO (UNSTRESSED) 

INITIAL CONDITIONS 



TOTAL NUMBER OF EQUILIBRIUM EQUATIONS 

, 

14 



NUMBER OF NON-ZERO STIFFNESS TERMS 

- 

105 



NUMBER OF EIGENVALUES BELOW SHIFT 

- 

2 



LINEAR STATIC CASES 



01:48:16 


1 USING STIFFNESS AT ZERO (UNSTRESSED) INITIAL CONDITIONS 



TOTAL NUMBER OF CASES TO SOLVE 

_ 

1 



NUMBER OF CASES TO SOLVE PER BLOCK 


1 



LINEAR STATIC CASES TO BE SOLVED: 





CASE: DEAD 





{ANALYSIS COMPLETE 



2007/05/03 01:48:16 

H. 


Ok > Aparece inmediatamente 
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6.- REVISAR LOS RESULTADOS 

Show Forces/Stresses > Frames/Cables 


SAP2000 vl 0.0.1 Advanced - tutorial. 2 



File Edit View Define Bridge Draw Select Assign Analyze Disglay Design Options Help 


l|OcS|H»| « o.| f 

| | a „i .| *h H 12 I •« S I i-afr. 


JS> .© ,© ,0 ,© I B 3-d » a js nu ♦ ?,J0 V 





Right Click on any joint for displacement values 


Start Animation I <^> I ■={> |GL0BAL ▼ II Kip. in, F - 
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Component > Axial Force [fuerza axial] 

Options > Show Values on Diagram [mostrar valores en el diagrama] 


Ok 



1 11 | 


SAP2000 vl 0.0.1 Advanced - tutorial 3. 2 


File Edit View Define Bridge Draw Select Assign Analyze Display Design Options Help 


||Da;|y^| « -■ f \ 


. nhH- . 


• & ,© & ,© 




£ S* H Itl 




'M Axial Force Diagram (DEAD) 


QEiS SS 





Right Click on any Frame Element for detailed diagram 3= | <=J> | GLOBAL ▼] Kip in, F •» 
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Los resultados aparecen en la ventana izquierda en la vista 3d, para observarlos mejor 
cambiaremos a la vista xz dando clic. 


K 


BBBI 


SAP2000 vl 0.0.1 Advanced - tutorial 3. 2 


File Edit View Define 

D B Stf 


Bridge Draw 

f a 

II Si SH i 


Select Assign Analyze Display Design Options Help 

p p p p p H td » .i n o- ac 

i | ^ | i. Il 8™ Set XZ )/cw I 




I 


. n o H - JJ 


'M Axial Force Diagram (DEAD) 


BBS X ■ -Z Plane 


Jin'S 



BBB1 


SAP2000 vl 0.0.1 Advanced - tutorial 3. 2 


File Edit View Define Bridge Draw Select Assign Analyze Display Design Options Help 


.| |1 n ff H~| .| 


D E? Q!tf 


I ► | p P P p P © | HI | 3-d »* b )! r» a AC I O O I ?n 0 I ■ 


'M Axial Force Diagram (DEAD) 


BBS 



% X-Z P 





Right Click on any Frame Element for detailed diagram 


<H*I GLOBAL - Kip, in. F - 
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Como podemos observar los resultados son iguales o muy parecidos a los realizados 
manualmente. 


Manualmente 


En el Sap 


18.0 kips 

24.0 kips 
14.40 kips 


18.67 kips 
24.82 kips 
14.74 kips 


3.2 ELEMENTO FINITO VIGA 

El elemento finito viga tiene eje recto con dos nodos y longitud L, existen dos grados de 
libertad por nodo: los desplazamientos y giros nodales. 



L 

x-0 x -L 

Fig. 3.9 Desplazamientos y giros nodales. 


Considerando que v, es el desplazamiento en el nodo “i” normal al eje longitudinal de la 
barra en el sentido y; es el giro del eje longitudinal de la barra en el nodo “i”. Definido 
como: 



•SI’S 

II 

(3.45) 

por definicion: 



Momenta: 

m z - { -)eP: 

(3.46) 

Cortante: 


(3.47) 

Carga: 

fcl 

II 

(3.48) 

P 

Para un elemento sin carga: — 

El 

= d ' v ,=o 

dx 4 

(3.49) 

Asi, la ecuacion diferencial que gobiema el comportamiento de una viga a flexion es: 




(3.50) 


Se propone como solucion a la ecuacion diferencial, cuya incognita es el desplazamiento, 
un polinomio cubico, tal que: 

v* = a x + a 2 x + a 3 x 2 + a pc’ (3.51) 


derivando, 
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dv 

dx 


= a 2 + 2a 3 x + 3a 4 x~ 


(3.52) 


se calcula la derivada del desplazamiento ya que esta relaciona a v con la otra incognita 

nodal: 0 = ^—. 

dx 

En las ecuaciones (3.51) y (3.52), las incognitas son ai, ci 2 , a j y a .4 (4 incognitas) por lo que 
se necesitan dos ecuaciones mas para determinarlas; los valores de ai, ci 2 , a , 3 y a .4 se 
obtienen de las condiciones de frontera, ya que al sustituir dos ecuaciones en dos nodos, 
tendremos cuatro ecuaciones y cuatro incognitas. 

Evaluando (3.51) en las condiciones de frontera, 
en x = 0 v* = a { 

x = L v* = + a 1 L + a 3 L 2 + a 4 L 3 

y evaluando (3.52) en las condiciones de frontera, 
dv * 


x = 0 
x = L 


dx 
dv * 
dx 


= a. 


— ci 2 + 2a 3 Z + 3 ci 4 L~ 


Ademas en los apoyos los desplazamientos deben ser igual a: 


y el giro nodal a: 


x = 0 
x = L 

x = 0 


v = v, 


V =v n 


dv 

dx 


= 0 , 


x = L 

Igualando ambas condiciones tenemos: 


ch L = n 

dx 


En forma matricial: 



E 

= a, 






V 2 

= a. 

+ &2 L 4 

• 4 

a A U 


0 

= a 2 






0 2 

= a~ 

, + 

2 cuL 

< 4 3a 4 L 

Vi 


“1 

0 

0 

0 " 

a ] 

*1 


0 

1 

0 

0 

a 2 

V 2 


1 

L 

L 2 

z 3 

a 3 



0 

1 

2 L 

3 L 2 _ 

_Cl 4 _ 

de 

libertad 

en 

los nodos. 


incognitas. 


a x — Vj 
a 2 = 


(3.53) 

(3.54) 
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Sustituyendo en: 
tenemos, 

Sustituyendo en: 
tenemos, 

Ordenando (3.55) y (3.56), 


v 2 — Ui + u,Z + a-^Is + a 4 L 

v, = v l + 6 X L + a 2 Is + a 4 L 3 
9 2 =a 2 + 2 a } L + 3 a 4 L 2 

0 2 = 9 X + 2 ci 4 L + 3a 4 L 
v 2 - Vj - 6 { L = « 3 Z 2 + « 4 Z 3 


(3.55) 


(3.56) 


9 2 -9\= 2 a 3 L + 3 a 4 L 2 

2 

Multiplicando la primera por y sumando tennino a tennino con la segunda, 

Z 

— — [v 2 -v, - 9 X L\ = -^/a 2 L - 2a 4 L 2 
L 

9 2 —9 x = 2 j/ 2 L + 3 a 4 L 2 

La suma es: 

— — [v 2 — v t - 9 X L\ + 9 2 — 9 X = a 4 L 2 
Z 


Despejando a 4 - 

a 4 




% 

L 2 


Finalmente, 


« 4 


2v 2 2v, 26> 9 2 
L 3 L 3 Is L 2 


(h 

L 2 


2v 2 2vi 0.0, 
Z 3 Z 3 Is Is 


(3.57) 


Sustituyendo (3.57) en: 


despejando a 3 , 


n 3 


^2 

2Z 2Z 


9 2 -9 x = 2 a 3 L + 3a 4 L 2 


3 

2 


« 4 Z 


n 3 


^__^L_3Z _2v^ 2^ ^ ^ 

zl ir + ir + i 2 + i 2 


a = ^i , 3v 2 3v i 3< ?i 3fl 2 

3 2Z 2Z Z 2 Z 2 2Z 2Z 

Finalmente, 
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<J 3 = -2 


L 


9 2 3vj 3 v 2 
~L Z r+ Z F 


(3.58) 


Sustituyendo en la aproximacion a la solution los valores de (3.53), (3.54), (3.57) y (3.58), 
en: 

v* = «[ + a 2 x + a 3 x 2 + a 4 x 2 


tenemos: 


v = v, + 9 x x + 


26*j 0 2 3v 1 3v, 

~L ~L 7F + 1F 


x 2 + 


0 t 0 2 2v, 2 v 2 

i 2+ i 2+ ir ir 


X 


Tomando tenninos comunes: 


6 bien: 
Donde: 


* 

, 3x 2 

2x 3 " 

+ 6 X 

2x 2 x 3 ’ 


_ 3x 2 

2x 3 ” 

+ 9 2 

2 3 

X X 

V = V, 

1 

+ , 

x + — 

+ v 2 



+ — 

1 

Z 2 

Z 3 _ 

i 

Z L 1 _ 

z 

Z 2 

Z 3 

z 

1 

N 

1 


v* = v l N l + 9 { N 2 + v 2 7V 3 + 9 2 N 4 


" 3x 2 2x 3 " 

1 — + — — 

; n, = 

_ 3x 2 2x 3 " 

Z 2 Z 3 

* 3 

Z 2 Z 3 


(3.59) 


2x 2 

x 3 

; n 4 = 

X 2 

X 3 

x 

H 

— 

-1 7- 

Z 

z 2 


Z 

Z 2 


Son las funciones de forma o de interpolation, tambien son conocidas como polinomios de 
Hennite. 

Las funciones de fonna o de interpolacion tienen las propiedades de valer uno al ser 
evaluadas en su nodo y valen cero en cualquier otro nodo. Asi para el grado de libertad 1, 


Aj (jc = 0) = 1 

N l (x = L) = 0 


para el grado de libertad 2, 

dNj . 4x 3x 2 

— - = l + — r 

dx L L 


para el grado de libertad 3, 


dN 2 

dx 

dN 2 

dx 


(x = 0) = 1 

(x = Z) = 0 


A 3 (x = 0) = 0 


N 3 (x = Z) = 1 

para el grado de libertad 4, 

dN 4 _ 2x 3x 2 
dx L L 2 
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dN 4 

dx 

dN 4 

dx 


(x = 0) = 0 

{x = L) = 1 


A continuacion se muestra la variation de las funciones de forma en el elemento 



0 0.5L L 


Fig. 3.10 Variation de las funciones de forma en el elemento. 


3.2.1 MATRIZ DE RIGIDEZ ELEMENTAL, ELEMENTO VIGA 



Fig. 3.11 Giros y desplazamientos nodales elemento viga. 

La ecuacion diferencial que gobierna el comportamiento a flexion del elemento viga es: 

r 4 

EI—j = 0 0 <x<L (3.60) 

dx 

La solucion de la ecuacion diferencial se evaluara por medio del metodo de residuos 
pesados. 

1) Se propone una aproximacion a la solution de la forma: 

4 

v* = Ajv, + N 2 0 { + A 3 v 2 + N 4 0 4 = Y, N i a i (3.61) 
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2) Sustituyendo la aproximacion en a la ecuacion diferencial: 


3) Se pesa el error o residuo, 


d A v 

El — - = R 

dx 4 


f R ■ WdO. = 0 

Jq ' 


(3.62) 


(3.63) 


\ L EI^\ ■W i dx = 0 

Jo /7r 4 


Integrando por partes, con: 


(3.64) 


u = W, 


1 4 * 

dv = — \-dx 
dx 4 


Sustituyendo, en (3.64), 


Elf ■ Wdx = EIW 

J ' rlE 


dx 0 


dx- 0 


(3.65) 


Para la segunda integral de (3.65), 


integrandola por partes, con: 

dW. 


El f EL 
l dx dx 


dv= d E- dx 
dx 


(3.66) 


Sustituyendo e integrando, 


r f dV [ dW, 1 , rT \dW t ~\ d 2 v L 
— T — ~ dx = EI — - — — 
i dx dx dx dx" 0 


l 

-El) 

0 J 


dv d 2 W : 


dx dx 


(3.67) 


Sustituyendo (3.67) en (3.65), 


d 3 v* L f dW - 1 J ? v* L r 

E/W — 5- -El ^ — - +El\ 

dx 3 0 dx dx" 0 J 




dx 


■r\dx = 0 (3.68) 


Despejando la integral definida, 


ItdVYd'w. 


dW, 1 d 2 v \\ L 


r <3 v u vv t avv i a v " ... a v - 

E7 — r 3 - dx = El — — -El W; — r- (3.69) 

J dx" dx 2 ^ ^ 2 0 ' //v 3 0 v 


dx \\ dx 0 


dx P 
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Como: 


= tfjVj + NA + N 3 v 2 + N 4 0 4 = £ N i a i 


dx 1 ~ dx 2 [ h 

^[XnA^^Nm+na+n^+na] 


(3.70) 


^!£ N A = ^ Vi + ^ 0 i+ ^ Vi+ ^± 02 

dx 1 [, =1 j dx dx" dx dx 


que se puede escribir, 


d 2 ff ^ lf jX d 2 N 2 d 2 N 3 d 2 N 4 1 6> 
dx 2 I “f ' ' J dx 2 dx 2 dx 2 dx 1 


(3.71) 


Por el Metodo de Galerkin, las fonciones de peso son las funciones de forma, 

W i = N i 
dW i _ dN . 
dx dx 
d 2 W t _ d 2 N t 
dx 1 dx 2 

Sustituyendo en (3.69), 


Md 2 N,J d 2 


dx" dx 


EE X 


g? 2 v d 


dx" I dx 0 


Sustituyendo (3.71), y variando i=l, 


E , ^ h u = 4^i®L'_4rfv L/ 

“ [_ dx J [_ dx dx dx" dx J v, |_ dx j dx 0 |_ <7x J 0 


para i = 2, 


Md 2 N,Jd 2 N, d 2 N , d 2 N-, d 2 N 4 ~]\0, , „Jrf 2 v*"| l £ _["rfV 


J 0 [_ dx J [_ dx dx dx dx" J v 2 


dx = El 


dx I dx 1° I dx 


para i = 3 
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L 

£/J 

0 


d 2 N 3 
dx 2 


d~N l d~N 2 d~N 3 d~N 4 


dx 2 dx dx dx 2 


para i = 4 


L, 

EI \ 


o L 


d~N A 

dx 2 


d N, d d Nt d~N, 


dx 2 dx dx dx 2 


4 


0, 


a 


e. 


dx = El 


1 

* 

> ( 

1 

dN 3 

L 

i 

* 

m 

1 

n 3 

3 

-El 

|_ dx 

dx 

0 

dx 



dx = El 


[-72 *1 

a v 

dN 4 

L 

i 

> 

* 

i 

n 4 


-El 

dx 2 

dx 

0 

dx 



Agrupando terminos: 


d 2 N x 
dx 2 
d 2 N , 


dx 2 

d 2 N 4 

dx 2 


r dx' 

El 

’ d 2 N, 

d 2 N 2 

d 2 N 3 d 2 N 4 " 

0 d 2 N 3 


dx 2 

dx 1 

dx 2 dx ' 


v i 


>i ' 

< 

> dx = 

M\ 

V 2 


v 2 

A 


m 2 


que tiene la forma: 

V 

en donde la matriz de rigidez elemental es: 

L 

[F] = J[5] 7 [D][B]Jx 


(3.72) 


(3.73) 


L, 

EI \ 


d 2 N l d 2 N x 

d 2 N x d 2 N 2 

d 2 N { d 2 N 3 

d 1 N l d 2 N 4 



dx 2 dx 2 

dx 2 dx 2 

dx 2 dx 2 

dx 1 dx 2 




d 2 N 2 d 2 N 2 

d 2 N 2 d 2 N 3 

d 2 N 2 d 2 N 4 






dx 2 dx 2 

dx 2 dx 2 

dx 2 dx 2 

dx 




My 

Simetrica 


d 2 N 3 d 2 N 3 

d 2 N 3 d 2 N 4 




v 2 


dx 2 dx 2 

dx 2 dx 2 


e 2 


m 2 


d N 4 d N 4 
dx 2 dx 2 


(3.74) 


Las integrales se resuelven anallticamente o numericamente por medio de la cuadratura de 
Newton-Cotes. 

EJEMPLO 3.3 
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Evaluar el termino k n de la matriz de rigidez elemental (3.74). Se necesita, 

3 ' 




, 3x 2 2x i 

1 ^ + 


L 1 


Is 


y sus derivadas: 

dN. 6x 6x 2 

L = 7 + — r 

dx Is L 

d 2 N { _ 6 12x 

~dx r ~~i 2+ ir 

Entonces: 


k 


ii 



d 2 N x d 2 N l 
dx 2 dx 2 




1 2x 

IF 


6 12x 

Is + Z 3 


dx 



L j 

EI\ 


o L 


36 144x 144x 
Z 4 Z 5 + ‘ 


Z 6 


dx 


= EI\ 


36 72x 2 48x 3 


— X 7 

z 4 z 5 


- + - 


z 6 


-l| o 


= EI 


36 (z) 72(Z) 2 | 48(Z) 3 


Z 


Z 3 


z° 


= El 


k 


1 1 


36 _ 72 48 

z 3 z 3 + z 3 

\2EI 

~Y~ 


EJEMPLO 3.4 

Evaluar el termino /c l2 de la matriz de rigidez elemental (3.74). 


L, 

k u = Z7 J 


d 2 N x 

dx 2 


d 2 N 2 

dx 2 


L, 

dx = Zzj 


_ 6 _ 12x 
Z 2 + Z 3 


4 6x 
Z + Is 


dx 


L 

EI \ 

0 


24 36x 48x 72x 2 

if IF 1F + ^lF 


ax 


102 



FUNDAMENTOS BASICOS DEL ELEMENTO F1NITO Hernandez Barrios Hugo FESA-UNAM 


= EI 


24a 18a 2 24a 2 24a 3 




L 4 L 4 + L 5 


= EI 


24a 42a 2 24a 3 


1} 


L 4 


- + - 


U 


= EI 


24(Z) 42(Z) 2 24(Z) 3 


Li 


V 


Is 


= El 


24 _42 24 

¥ ¥ + ¥ 


, 6EI 
kn= ~¥ 


Los demas terminos se obtienen de forma similar, finalmente la matriz elemental de rigidez 
es, 


(3.75) 



" 12 

6L 

-12 

6L ' 

v i 

El 

6L 

4 L 2 

-6 L 

2li 

4 

¥ 

-12 

-6L 

12 

-6L 

V 2 


6L 

2 Is 

-6 L 

4 Is 

A. 


EJEMPLO 3.5 

Determinar los elementos mecanicos en la siguiente viga. 

y 

A 

500 


2000 


i 


■>X 


10 ’ , 10 ’ 

l< >l< >\<r 


12 ’ 


->1 


Fig. 3.12 Viga del ejemplo 3.5. 


2000 


500 


n r* 

1 1 ' 2 

3 

V 

4 

j 10’ 10’ 

& 

12’ 


k >< 

-> j < 



Fig. 3.13 Topologia del modelo. 


ELEMENTO 

CONECTIVIDAD 

MODULO DE 

MOMENTO 
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ELASTICIDAD 

DE INERCIA 

1 

1-2 

E 

I 

2 

2-3 

E 

I 

3 

3-4 

E 

I 


Las matrices de rigidez elemental, segun (3.75) son: 

ELEMENTO 1 



L 

% 

V 2 

0 2 



" 12 

60 

-12 

60 " 

v i 

M= EI 

L J 1000 

60 

400 

-60 

200 

0, 

-12 

-60 

12 

-60 

V 2 


60 

200 

-60 

400 

0 2 

ELEMENTO 2 






V 2 

0 2 

V 3 




' 12 

60 

-12 

60 ' 

V 2 

[K]= EI 

L J 1000 

60 

-12 

400 

-60 

-60 

12 

200 

-60 

0 2 

V 3 


60 

200 

-60 

400 

^3 


ELEMENTO 3 



^3 

0, 

v 4 

0 \ 


6.94 

41.7 

-6.94 

41.7 " 

V 3 

41.7 

333.0 

-41.7 

167 

02 

-6.94 

-41.7 

6.94 

-41.7 

^4 

41.7 

167 

-41.7 

333.0 

0, 


La matriz de rigidez total, referida a un sistema global es, 


a 


0 , 


0 ' 


0 , , 
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12 

60 

-12 

60 





Vi 


60 

400 

-60 

200 





4 


-12 

-60 

12 + 12 

-60 + 60 

-12 

60 




El 

60 

200 

-60 + 60 

400 + 400 

-60 

200 



e 2 

1000 



-12 

-60 

12 + 6.94 

-60 + 41.7 

-6.94 

41.7 

V 3 




60 

200 

-60 + 41.7 

400 + 333 

-41.7 

167 

0, 






-6.94 

-41.7 

6.94 

-41.7 

V 4 






41.7 

167 

-41.7 

333 

o* 


Se tiene entonces, 



' 12 

60 

-12 

60 


60 

400 

-60 

200 


-12 

-60 

24 

0 

El 

60 

200 

0 

800 

L 

0 

0 

-12 

-60 


0 

0 

60 

200 


0 

0 

0 

0 


0 

0 

0 

0 


0 

0 

0 

0 


v i 


V,' 

0 

0 

0 

0 


3 


M, 

-12 

60 

0 

0 


v 2 


C 

-60 

200 

0 

0 


d 2 

> = ' 

m 2 

18.9 

- 18.3 

- 6.94 

41.7 




V, 

- 18.3 

733 

- 41.7 

167 


£ 3 


m 2 

- 6.94 

- 41.7 

6.94 

- 41.7 


V 4 


V A 

41.7 

167 

- 41.7 

333 


.* 4 . 


m 4 


Con las siguientes condiciones de frontera, para el empotre, 

v>=0 

e x = o 

en el apoyo movil, v 3 = 0 

Condensando la matriz global, queda, 



"24 

0 

60 

0 

0 


V 


\V 2 ] 


[ v 2 1 


0 

800 

200 

0 

0 


e 2 


m 2 


-500 

El 












60 

200 

733 

-41.7 

167 


o, : 


M, 


M, 

1000 













0 

0 

-41.7 

6.94 

-41.7 






-2000 


0 

0 

167 

-41.7 

333 


kJ 


kJ 


. M, j 


Ya que, M 2 = -500 y V 4 = -2000 
Resolviendo el sistema de ecuaciones obtenemos, 
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y finalmente, 


' 

V 2 


' 304348 

e 2 

i 

29810 

o. 

> = — 

EI 

-121739 

V 4 


-2636991 

A. 


-269166 

v. 9 


' > 

Vi 


' 0.00000 ' 



0.00000 



304348 

e 2 

i 

29810 

> 


~ eA 

0.0000 

o. 


-121739 

V 4 


-2636991 

A, 


-269166 


Para el calculo de momentos y cortantes, de la relation (3.75): 



“ 12 

6L 

-12 

6L ' 

Vi 



EI 

6L 

4L 2 

-6L 

21} 

4 


M x 

1} 

-12 

-6L 

12 

-6L 

V 2 

— ■= 

v 2 


6L 

2 Is 

-6 L 

1 

A. 


M , 

L ^9 


ELEMENTO 1 



" 12 

60 

-12 

60 " 

0 


>i ' 

EI 

60 

400 

-60 

200 

0 


M x 

1000 

-12 

-60 

12 

-60 

304348/ EI 

’ — < 

v 2 


60 

200 

-60 

400 

29810 / EI 

v. 9 


m 2 

L z 9 


vA 


-1864' 

Mi 

> = < 

-12299 

> 

v 2 


1864 

m 2 _ 


-6337 

^ 9 


ELEMENTO 2 



" 12 

60 

-12 

60 " 

' 304348/ EI 


'vA 

EI 

60 

400 

-60 

200 

29810 /EI 


m 2 

1000 

-12 

-60 

12 

-60 

0.00 

> — ■= 

V, 


60 

200 

-60 

400 

-121739 / EI 


M, 

< J 9 


106 



FUNDAMENTOS BASICOS DEL ELEMENTO F1N1TO Hernandez Barrios Hugo FESA-UNAM 


V 2 ' 


-1864' 

m 2 

> = < 

5837 

> 

V 3 


1864 

m 3 _ 


- 24473 


ELEMENTO 3 



" 6.94 

41.7 

-6.94 

41.7 “ 

0.0 


’V 3 ' 

El 

41.7 

333.0 

-41.7 

167 

-121379/E/ 


m 3 

1000 

-6.94 

-41.7 

6.94 

-41.7 

-2636991/E/ 

=■ — ' 

v 4 


41.7 

167 

-41.7 

333.0 

-269166 /El 

v. 9 


m 4 


> 3 ' 


2000 ' 

m 3 

> = < 

24473 

v 4 


2000 

M 4. 


— i 66 


Para calcular las reacciones: 



“ 12 

60 

-12 

60 

0 

0 

0 

0 


0 


*1 


60 

400 

-60 

200 

0 

0 

0 

0 


0 


Mi 


-12 

-60 

24 

0 

-12 

60 

0 

0 


304348 


0 

El 

60 

200 

0 

800 

-60 

200 

0 

0 

1 j 

29810 

„ _ J 

-500 

L 

0 

0 

-12 

-60 

18.9 

-18.3 

-6.94 

41.7 

El 

0 

1 — 1 



0 

0 

60 

200 

-18.3 

733 

-41.7 

167 


-121739 


0 


0 

0 

0 

0 

-6.94 

-41.7 

6.94 

-41.7 


-2636991 


-2000 


0 

0 

0 

0 

41.7 

167 

-41.7 

333 _ 


-269166 

l j 


0 


de donde: 

R, =-1864 
Mj = -12299 
R 3 = 3864 

3.3 MATRIZ DE RIGIDEZ ELEMENTAL EN COORDENADAS GLOBALES 

Anteriomiente se ha encontrado la matriz de rigidez en coordenadas globales para un 
elemento viga con el eje recto, tal como se mostro en la expresion (3.75). 

Se considera que esta referida a coordenadas globales ya que los ejes locales y globales 
coinciden. La matriz de rigidez elemental que considera carga axial, cortante y momenta, se 
obtiene de sumar la matriz de rigidez del elemento annadura y la del elemento viga, 
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“l V 1 


EA 

L 

0 

0 

EA 

L 

0 

0 


0 

12 El 

ff 

6EI 

FF 

0 

12 El 

nr 

6EI 

IF 


% 

0 

6EI 

IF 

4EI 

L 

0 

6EI 

'IF 

2 El 
L 


u 2 

EA 

L 

0 

0 

EA 

L 

0 

0 


v 2 0 2 

0 0 

12 El 6 El 

L 3 FF 

_ 6EI 2 El 

L 2 L 

0 0 

12 El _6EI 
L 2 L 2 

_ 6EI 4 El 

L 2 L 


M, 

Vi 

Ox 

u 2 

V 2 

e 2 


(3.76) 


donde: u l3 v, , 0 ] , u 2 , v 2 y 0 2 son los grados de libertad del elemento referidos a un 

sistema coordenado local. Para el caso de que el eje de la viga este inclinado un angulo, 6, 
con respecto a la horizontal, la matriz de rotacion que pennite transfonnar la matriz de 
rigidez a un sistema global es: 


c s 0 0 0 0 

s —c 0 0 0 0 

0 0 1 0 0 0 

0 0 0 c s 0 

0 0 0 —s c 0 

0 0 0 0 0 1 


(3.77) 


que relaciona los grados de libertad en coordenadas locales u { , v, , 9 X , u 2 , v 2 y 0 2 con los 
referidos a coordenadas globales: u[, vj , 0[, u 2 , v 2 y 0 2 . 



Fig. 3.14 Direccion positiva de los ejes locales x-y; y los globales x -y . 
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u. 


c 

s 

0 

0 

0 

0 " 


Mj 

L 


s 

-c 

0 

0 

0 

0 


v l 

% 


0 

0 

1 

0 

0 

0 

< 

0 [ 

> 



0 

0 

0 

c 

s 

0 


U 2 

V 2 


0 

0 

0 

-s 

c 

0 


v 2 



0 

0 

0 

0 

0 

1 


& 


(3.78) 


De tal manera que la matriz de rigidez elemental referida a coordenadas globales esta 
definida por: 




es decir: 


c 2 e + s 2 f 

cs(e-f) 

~sg 

-e 2 e-s 2 f 

cs(f-e) 

~sg 

cs(e-f) 

s 2 e + c 2 f 

eg 

es(f-e) 

-s 2 e-c 2 f 

eg 

~sg 

eg 

h 

sg 

-eg 

h 

2 

-c 2 e-s 2 f 

cs(f-e ) 

sg 

c 2 e + s 2 f 

cs(e-f ) 

sg 

cs(f-e) 

-s 2 e-c 2 f 

-eg 

cs(e-f ) 

s 2 e + c 2 f 

-eg 

~sg 

eg 

h 

2 

sg 

-eg 

h 


donde: 

EA , 12 El 6EI , 4EI 

e = — , f = — — , g = — — , h = y s = send ; c = cos 6 

L L L L 


(3.79) 


Las fuerzas intemas en la barra en coordenadas locales se determina con: 
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EA 

L 

0 

0 

EA 

L 

0 

0 


0 

12 El 

~F~ 

6EI 

FF 

0 

12 El 

FT 

6EI 

IF 


0 

6EI 

IF 

4 El 
L 

0 

6EI 

'IF 

2 El 
L 


EA 

L 

0 

0 

EA 

L 

0 

0 


0 0 

12 El 6EI 

E L 2 
6EI 2 El 

FF FF 

o o 

\2EI _6EI 

L 2 L 2 

6EI 4 El 

FF FF 


M| 

V 1 

u 2 

V 2 

0 2 


(3.80) 


6 bien, una vez conocidos los grados de libertad en coordenadas globales, 


y i 


e\ 


y 2 


EA 

0 

0 

EA 

0 

0 



















L 



L 











0 

\2EI 

6EI 

0 

\2EI 

6 El 



0 

0 

0 

0" 


r i b 

L 2 

L 2 

L 2 

L 2 

c 

s 


M, 

0 

6EI 

4EI 

0 

6EI 

2EI 

s 

-c 

0 

0 

0 

0 


F 

L 2 

L 

L 2 

L 

0 

0 

1 

0 

0 

0 


o\ 

EA 

0 

0 

EA 

0 

0 

0 

0 

0 

c 

.V 

0 


Ur , 







L 



L 



0 

0 

0 

-s 

c 

0 


V 2 

0 

12 El 

6EI 

0 

\2EI 

6EI 

0 

0 

0 

0 

0 

1 



L 2 

L 2 

L 2 

Is 











0 

6EI 

2EI 

0 

6EI 

4EI 






















L 2 

L 


L 2 

L 










(3.81) 


Si existe una carga uniformemente distribuida en el miembro, como la de la figura (3.16), el 
vector de cargas nodales es: 


M=(-) M=(+) 



z 


Fig. 3.15 Sentido positivo de los ejes locales. 
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wL wL 



Fig. 3.16 Cargas nodales equivalentes. 


f F e } = 

0 ~'± 

o 

1 

wL 

wL 2 

( ) 

2 

12 

~Y 

12 


El vector de cargas referido a un sistema global es: 

[F‘] = [rf{F'} 

Los valores de [ J se agregan al vector de carga global; las cargas concentradas y los 
momentos concentrados se agregan al vector de carga nodal. 

Al agrupar las cargas y las rigideces, se obtiene el sistema de ecuaciones: 

(3.82) 

EJEMPLO 3.6 

Determinar los desplazamientos y rotaciones de los nodos del marco mostrado en la figura 
(3.17), considerando los siguientes datos: 
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w=10Q kg/m 



Fig. 3.17 Marco del ejemplo 3.6. 


0'. *■ I A 


e; 


m 


TXT 


o\- . 


m 


9L 


I 

1 4 

Fig. 3.18 Topologia del marco del ejemplo 3.6. 


1 ,400 kg 


L 


1 ,400 kg 


933.33 kg-m 


933.33 kg-m 


Fig. 3.19 Cargas nodales equivalentes en el elemento 2. 

Datos: 

Concreto clase 1, (NTC-RCDF2004) con un f ' c = 250-^- , por lo que 

cur 

E = 14000^250 = 221359^ = 221.359— ^ . 

cm cm 

Columnas de 50x50 cm , trabe de 30x40 cm . 


Elemento 

A (cm 2 ) 

I (cm 4 ) 

L (cm) 

1 

2500 

520834 

250 

2 

1200 

160000 

400 

3 

2500 

520834 

250 
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Conectividades: 


Constantes: 


Elemento 

nodo 

i 

.i 

1 

1 

2 

2 

2 

3 

3 

4 

3 


Elemento 

EA 
e = — 

L 

12 El 

Z 3 

6EI 

g= e 

, 4 El 

h = 

L 

c = cos# 

s = senO 

1 

2214.0 

88.54 

11068 

1844661 

0.0 

1.0 

2 

664.0 

6.64 

1328.15 

354174 

1.0 

0.0 

3 

2214.0 

88.54 

11068 

1844661 

0.0 

1.0 


Matriz de rigidez elemental: 


ELEMENTO 1 



ELEMENTO 2 



ELEMENTO 3 



u\ 

v i 


u' 2 

V 2 

4 


" 88.54 

0 

-11068 

-88.54 

0 

-11068 “ 

Mj 

0 

2214 

0 

0 

-2214 

0 

vi 

-11068 

0 

1844661 

11068 

0 

922330.5 

o[ 

-88.54 

0 

11068 

88.54 

0 

11068 

U 2 

0 

-2214 

0 

0 

2214 

0 

V 2 

-11068 

0 

922330.5 

11068 

0 

1844661 

0 2 

u 2 

V 2 

#; 

u 3 

V 3 

% 


664 

0 

0 

-664 

0 

0 

U 2 

0 

6.64 

1328.15 

0 

-6.64 

1328.15 

V 2 

0 1328.15 

354174 

0 

1328.15 

177087 

# 2 

-664 

0 

0 

664 

0 

0 

m 3 

0 

-6.64 

-1328.15 

0 

6.64 

-1328.15 

V 3 

0 1328.15 

177087 

0 

1328.15 

354174 


u 4 

v 4 

^4 

u 3 

V 3 



88.54 

0 

11 068 

-88.54 

0 

-11068 " 

« 4 

0 

2214 

0 

0 

-2214 

0 

^4 

-11068 

0 

1844661 

11068 

0 

922330.5 

0 \ 

-88.54 

0 

11068 

88.54 

0 

11068 

u 3 

0 

-2214 

0 

0 

2214 

0 

V 3 

-11068 

0 

922330.5 

11068 

0 

1844661 

0\ 
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Matriz de rigidez global: 


u \ 

v| 

4 ' 

u 2 

v 2 

e\ 


V 3 

4 

« 4 




" 88.54 

0 

-11068 

- 88.54 

0 

-11068 

0 

0 

0 

0 

0 

0 

«; 

0 

2214 

0 

0 

-2214 

0 

0 

0 

0 

0 

0 

0 


-11068 

0 

1844661 

11068 

0 

9223305 

0 

0 

0 

0 

0 

0 

4 

- 88.54 

0 

11068 

75254 

0 

11068 

-664 

0 

0 

0 

0 

0 

u 2 

0 

-2214 

0 

0 

222 Q 64 

132815 

0 

- 6.64 

132815 

0 

0 

0 

V 2 

-11068 

0 

9223305 

11068 

132815 

2198835 

0 

-132815 

177087 

0 

0 

0 

4 

0 

0 

0 

-664 

0 

0 

75254 

0 

11068 

- 88.54 

0 

11068 

m 3 

0 

0 

0 

0 

- 6.64 

-132815 

0 

222 Q 64 

-132815 

0 

-2214 

0 

4 

0 

0 

0 

0 

132815 

177087 

11068 

-132815 

2198835 

-11068 

0 

9223305 

4 

0 

0 

0 

0 

0 

0 

- 88.54 

0 

-11068 

88.54 

0 

-11068 

m 4 

0 

0 

0 

0 

0 

0 

0 

-2214 

0 

0 

2214 

0 

V 4 

0 

0 

0 

0 

0 

0 

11068 

0 

9223305 

-11068 

0 

1844661 

4 


La matriz global de rigideces condensada, considerando las restricciones de los apoyos es: 


u 2 

V 2 


u 3 

V 3 

^3 


752.54 

0 

11068 

-664 

0 

0 

tl 2 

0 

2220.64 

1328.15 

0 

-6.64 

1328.15 

V 2 

11068 

1328.15 

2198835 

0 

-1328.15 

177087 

d' 2 

-664 

0 

0 

752.54 

0 

11068 

« 3 

0 

-6.64 

-1328.15 

0 

2220.64 

-1328.15 

^3 

0 

1328.15 

177087 

11068 

-1328.15 

2198835 

0 ; 



El sistema de ecuaciones lineales a resolver es: 

donde: 

{f 2 }=[ 2 -1.4 -93.33 0 -1.4 93. 33] 7 
las fuerzas en t y momentos en t-cm. 


“752.54 

0 

11068 

-664 

0 

0 

u 2 


+ 2 

0 

2220.64 

1328.15 

0 

-6.64 

1328.15 

V 2 


-1.4 

11068 

1328.15 

2198835 

0 

-1328.15 

177087 

0 : 

« 3 


-93.33 

-664 

0 

0 

752.54 

0 

11068 

> < 

0 

0 

-6.64 

-1328.15 

0 

2220.64 

-1328.15 

V 3 


-1.4 

0 

1328.15 

177087 

11068 

-1328.15 

2198835 _ 

03. 


93.33 

V J 


de donde se obtienen los valores nodales: 
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' , ' 

u 2 


' 0.0282 ' 

V 2 


-0.0004 

0' 2 

> — < 

-0.0001 

u 3 


0.0259 

V 3 


-0.0007 

A. 


-0.0000 

C J 


donde los desplazamientos estan en cm y los giros en radianes. 


' 88.54 

0 

-11068 

- 88.54 

0 

-11068 

0 

0 

0 

0 

0 

0 

0 

0 

2214 

0 

0 

-2214 

0 

0 

0 

0 

0 

0 

0 

0 

-11068 

0 

1844661 

11068 

0 

9223305 

0 

0 

0 

0 

0 

0 

0 

- 88.54 

0 

11068 

75254 

0 

11068 

-664 

0 

0 

0 

0 

0 

0.0282 

0 

-2214 

0 

0 

222 Q 64 

132815 

0 

- 6.64 

132815 

0 

0 

0 

- 0.0004 

-11068 

0 

9223305 

11068 

132815 

2198835 

0 

-132815 

177087 

0 

0 

0 

- 0.0001 

0 

0 

0 

-664 

0 

0 

75254 

0 

11068 

- 88.54 

0 

11068 

0.0259 

0 

0 

0 

0 

- 6.64 

-132815 

0 

222 Q 64 

-132815 

0 

-2214 

0 

- 0.0007 

0 

0 

0 

0 

132815 

177087 

11068 

-132815 

2198835 

-11068 

0 

9223305 

- 0.0000 

0 

0 

0 

0 

0 

0 

- 88.54 

0 

-11068 

88.54 

0 

-11068 

0 

0 

0 

0 

0 

0 

0 

0 

-2214 

0 

0 

2214 

0 

0 

0 

0 

0 

0 

0 

0 

11068 

0 

9223305 

-11068 

0 

1844661 

0 


Mi 

2 

-1.4 

-93.33 

0 

-1.4 

93.33 

K 

K 

_ m 4 

A continuacion se explica paso a paso como realizar el ejemplo del marco en el programa 
SAP 2000 V. 10. 

MARCO DEL EJEMPLO 3.6 [P.112.] 

1.- CREACION DE UN ARCHIVO NUEVO 

File > New Model 
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Initialize Model from Defaults with Units [seleccionar unidades] 
Kgf, m, C [kg fuerza, metros, centigradosl 



Dar clic en 2D Frames [marcos] 
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Number of Stories 1 [entrepisos] 

Number of Bays 1 [claros] 

Story Height 2.5 [m] [altura de entrepiso] 
Bay Width 4 [m] > Ok [ancho del clarol 



Inmediatamente aparece 
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Escoger un lado para trabajar, en este ejemplo trabajaremos del lado derecho [dar clic 
ventana derecha], 

2.- DEFINIR MATERIALES Define > Materials 



en 


118 





FUNDAMENTOS BASICOS DEL ELEMENTO F1NITO Hernandez Barrios Hugo FESA-UNAM 



Cone > Modify/Show Material 



□taal 


Display Color 


Material Name 


Type of Material 
(• Isotropic 


Type of Design 


C Orthotropic 
C Uniaxial 


Analysis Property Data 
Mass per unit Volume 
Weight per unit Volume 
Modulus of Elasticity 
Poisson's Ratio 


Design Property Data (ACI 31 8-05/1 BC 2003) 
Specified Cone Comp Strength, f'c 1 281 2278.5 
B ending R einf. Yield S tress, fy 1 421 841 78 

Shear Reinf. Yield Stress, fys |281 22785 


|2.531E+09 


I - Lightweight Concrete 


Coeff of Thermal Expansion 
Shear Modulus 


Advanced Material Property Data 


Time Dependent Properties. 


Cancel 


Material Damping Properties. 


Stress-Strain Curve Definitions. 


* * »" »- »» ^ ^ ^ ■»- ■» T , I 

B6g|H»|» ? S > J8 J9 p ^ J3 | BB| m w « w •» & \ o- » | °b a | ^ . | || I -[ ■ - 1 ,| j]n rr H-| ,| 


XZPI»>e@Y-0 X0.G0 YO.OO ZO.M 3 

Modulus of Elasticity >2 213 590 000 [Kg/ 111 2 ] 

Specified Cone Comp Strength, f c > 2 500 000 [Kg/m 2 ] > Ok > Ok 


Material Property Data 
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3.- DEFINIR SECCIONES 

Define > Frame Sections 
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Choose Property Type to Add > Add Rectangular 



Click to > Add New Property 
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Section Name > Col 50x50 [cm] [nombre de la seccion] 
Material > Cone [material] 

Depth > 0.5 [m] , Width > 0.5 [m] > Ok [profimdo, ancho] 
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% AP200C v 1 0 . 0 ■ Advanced - Untitled] 


File Edit View Define Bridge Draw Select Assign Analyze Dismay Design Options Help 


D BS H } & ’ P @ ,© 1® P H 3-d »). *z yz frf « O 0 X- „ 

I • a . nfr«- . 

J - i. . 



X-Z Plane @Y-0 



X0.23 YO.OO Z4 11 


■ GLOBAL _fl|K 




Choose Property Type to Add > Add Rectangular 


Click to > Add New Property 


% SAP200C vio.o • Advanced Untitled] 

BBS 

File Edit View Define Bridge Draw Select Assign Analyze Dismay Design Options Help 


Dtf Hii i t & > P i*> .© & P H M « n Oise * » fju ? . 

I • a . nhH- . 

|| » r„ ! 1 it A. .[ II c-i Hi If| 1 I i. ^ J 




PS 



X-Z Plane @Y=0 


~3 


X0.23 YO.OO Z4.11 


Section Name > Trabe 30x40 [cm] [nombre de la seccion] 


Material > Cone [material] 

Depth > 0.4 [m] , Width > 0.3 [m] > Ok > Ok [profundo, ancho] 
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4.- ASIGNAR CONDICIONES DE APOYO, TIPOS DE SECCIONES Y CARGAS 
4.1 Asignar Condiciones de Apoyo 
Seleccionar el primer apoyo con el mouse 
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Assign > Joint > Restraints 
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Apoyo empotrado > Dar clic > Ok 
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Seleccionar el segundo apoyo con el mouse 



Assign > Joint > Restraints 
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Apoyo empotrado > Dar clic > Ok 



4.2.- Asignar Tipos de Secciones 
Seleccionar las columnas con el mouse 
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Assign > Frame/Cable/Tendon > Frame Sections 
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Seleccionar Col 50x50 > Ok 



Seleccionar la trabe con el mouse 
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Assign > Frame/Cable/Tendon > Frame Sections 


% 



.0.1 Advanced (Untitled) 


□SB 


Edit 

View 

Define Bridge Draw Select Assign Analyze Disglay Design Options Help 


1 

□ & 


■o t S * @ & & ,© P H 3-d »y yz ; Ct«A O 0 Jf, . 

I - I 

□ - . nriH- . 

II : 


l \& 

- Mm? ** 



¥ 



X Frame Sections 





Global J|- . J 


X0.00 Y0.00 Z0.00 


Seleccionar Trabe 30x40 > Ok 
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5.3.- Asignar Cargas 

Seleccionar la trabe donde sera aplicada la carga uniforme 



Assign > Frame/Cable/Tendon Loads > Distributed 
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SAP2000 vl 0.0.1 Advanced -(Untitled) 


BBBI 


□ e 

Edit View Define 

Bridge Draw Select 

Assign Analyze Display 

Design Options Help 

dcs ay « 

t s * sb>s 

Frame/Cable/T endon 


xz yz nv Cf at | <► o | ?« 0 | 3& , | J|l*|D"|,|]jnr7M’ - |-| 

II S3 

• , _ 1 j < ; 

a «h y iii 



?| 

% D 

SdM 

> 

d]B Frame Sections 





| GLOBAL £]|Kgf. m.C 


X0.00 Y0.00 Z0.00 



Unifomi Load >700 [kg] [carga uniforme] 
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Seleccionar ahora el nodo donde sera aplicada la fuerza 


* 

d 

o 

> 

o 

p 

© 

Cl 

< 

Advanced - (Untitled) 

BBJB 

File 

Edit View Define Bridge Draw Select Assign Analyze Display Design Options Help 


Da: H a -i 

f @ m *> >* k o * 0 y a . x- 

U - . niiW- . 


5 si -i . 

j] L C- y Hi A - - x | i- S£ |i . 


IS 

% D 

- n]Q I 'M Frame Distributed Loads (DEAD) 

- nisi 



X-1.25YO.OOZ2.13 | GLOBAL _-J | Kgf. m. C _£] 


Assign > Joint Loads > Forces 
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Force Global X > 2000 > Ok [fuerza en el eje X global] 
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6.- ANALIZAR [EJECUTAR] EL PROGRAMA 

Analize > Run Analisis 
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Case Name > Modal > Click > Run/Do Not Run Case > Click [tipos de analisis] 
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Run Now > Click 

Save Model File As [crear un folder, poner nombre y guardarlo ahi] 
Save 
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Google 

Deskt... 


Google Tropa 

Deskt... Jambor... 


ffi 

Tropa 

ambor... 


Relaoon 

play* 


SM 0 

Adobe Reader RecordNow! 

7.0 

0 * 

Norton My Network 
SystemWorks Places 

I* & 

Tunes My Computer 

m a 

e-Wortertouc... AutoCAD 2004 




■Mm)i x 










LINEAR EQUATION SOLUTION 

FORMING STIFFNESS AT ZERO (UNSTRESSED) INITIAL CONDITIONS 

TOTAL NUMBER OF EQUILIBRIUM EQUATIONS = 12 

NUMBER OF NON-ZERO STIFFNESS TERMS - 7B 

NUMBER OF EIGENVALUES BELOW SHIFT = 0 

LINEAR STATIC CASES 

USING STIFFNESS AT ZERO (UNSTRESSED) INITIAL CONDITIONS 

TOTAL NUMBER OF CASES TO SOLVE - 1 

NUMBER OF CASES TO SOLVE PER BLOCK - 1 

LINEAR STATIC CASES TO BE SOLVED: 

CASE: DEAD 


Unused 

Tesis 





Deskto... 



= 1 

- 

0 


Papdera de 
reddaje pr... 


Ok > Aparece Inmediatamente 


K 


UsM I 


SAP2000 vl 0.0.1 Advanced - tutorial 3. 6 


File Edit View Define Bridge Draw Select Assign Analyze Display Design Options Help 


D & S & ' | 9 l[fl[ ► \& & & & & vy xz yz rw 

j a f w t | « .| !| it H 'dl-v | “"a | ±. an .| 


j a y a . 


I Q • 


nnH’ - 


| M Deformed Shape (DEAD) 
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7.- REVISAR LOS RESULTADOS 

Show Forces/Stresses > Joints 


SAP2000 vl 0.0.1 Advanced - tutorial 3. 6 


File Edit View Define Bridge Draw Select Assign Analyze Display Design Options Help 


1 a l, % | a „i ,| || a a nt y | ■» 2 . | i- <a & .1 
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Quitar la paloma en Show as Arrows > Ok 




3-D View 


| => || GLOBAL 


▼ 1 1 Kgf, m, C ▼ 
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Los resultados aparecen en la ventana izquierda en la vista 3d, para observarlos mejor 
cambiaremos a la vista xz dando clic. 



1 SHj 


SAP2000 vl 0.0.1 Advanced - tutorial 3. 6 


File Edit View Define Bridge Draw Select Assign Analyze Display Design Options Help 


□ n r? M • ij 


DS 0‘J 


H m w « n » 34f « o ?h 0 y a . 




Joint Reactions (DEAD) 








O I => (global ▼] | Kgf, m, C 
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Asi obtenemos los resultados de las fuerzas y momentos de los nodos 1 y 4 del ejemplo 3.6 
Nodo 1 Nodo 4 


FI =R1X = -337.66 kg 
F3 =R1Y = 3149.15 kg 
M2 = Ml =-1367.52 kg-m 


FI =R4X = -1662.34 kg 
F3 = R4Y = 3807.64 kg 
M2 = M4 = -2315.48 kg-m 


3.4 ELEMENTO TRIANGULAR LINEAL 

Un elemento triangular lineal tiene lados rectos y tres nodos, uno en cada esquina. El nodo i 
se ubica arbitrariamente pero debe seguirse una numeracion en el sentido contrario al 
movimiento del reloj . 

Los valores nodales de (j) son ® ,® ; y ® /t y las coordenadas nodales son: 

(w) 

*j,yj) 

Xk’Vk) 

Fig. 3.20 Coordenadas nodales. 



Fig. 3.21 Valores y coordenadas nodales. 

El polinomio de interpolacion es: 

(/) = a { + a 2 x + a 3 y (3.83) 

Las rigideces nodales son: 
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4> = < & t en a = x ; y = y t 

(f> = en x = x j y = y } 

<t> = en x = x k y = y k 

Sustituyendo: 

=a l + a 2 x i + a,y i 
0 ; = a } + a 2 Xj + a 3 y j 
O k =a 1 +a 2 x k +a 3 y k 

Resolviendo el sistema de ecuaciones simultaneamente: 


1 r 


2 A e L 


( x } y k ~ Wj ) °, + {x k y, - x ,y t ) + (wj ~ x jyj ) °* 


2 A e 


(y, -y k )®i+{y k - y t ) + [y t - y k ) ^a 


a, = 


= 4 l X k ~ X j )® , + (*/ - X k ) (I) J + (Xj ~ X i )® k ] 


2 A 


Donde el elemento, A es el area del triangulo. 


2 A e = 


1 x, y t 

1 *,• y t 

1 x u y k 


Sustituyendo en <j) = a x + a 2 x + a 3 y (3-83), 


flj 


a 2 


a 3 


1 

~iF 

i 

~2A e 

1 

~iF 


[«/®/ +a j®j +« a ° a ] 

-[A®, + 4®, + A®*] 

\rp,+rpj +/ a ° a ] 


Donde: 


(3.84) 

(3.85) 

(3.86) 

(3.87) 

(3.88) 

(3.89) 

(3.90) 

(3.91) 

(3.92) 


(3.93) 

(3.94) 

(3.95) 
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-ts: 

X 

1 

x'’ 

II 

(3.96) 

a i = ( x / c y, ~ x ,y k ) 

(3.97) 

$ 

?r* 

II 

"x" 

i 

^x 

(3.98) 

> 

ll 

i 

(3.99) 

?s 

I 

II 

(3.100) 

> 

ll 

yP 

i 

(3.101) 

Yi=( x k- x j) 

(3.102) 

ii 

yT 

i 

(3.103) 

II 

i 

_x: 

(3.104) 


Se puede observar que las funcioncs de forma son constantes una vez que se conocen las 
coord enadas nodales. 

En (j) = a x + a 2 x + a 3 y lo que significa que los valores nodales son funcion de a y y. 

+a j®j +ar *®*] + ^7[A®,- + Pj®j + Pfii + y^j +r k ® k ]y 

(3.105) 

Ordenando: 

<!> = 7^7 [«,- + Pt* + r,y] + 7^7 [ a j + Pj x + y,y] + 7^7 K + P k x + ® * 

(3.106) 

o bien: 

^ = 7V.® ; + Nj<t> j + N k <& k (3.107) 

Donde: 

N t = 7^7 [«/ + A* + (3 • 1 08) 

^ = 2Z ^ + PjX + (3 ' 1 ° 9) 

N k = 7^7 K + v] (3.110) 

Polinomios lineales una vez que se conocen a , [5 y y , con funciones de forma en funcion 
de x y y. 

Los gradientes — y — son constantes en todo el elemento, por ejemplo: 

dx dy 
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Pero: 


d(j) 

_dN t 

®, + 

dx 

dx 


dN . 

d 

1 

dx 

dx 

_2A e 

dN t 

_ Pi 


dx 

2 A 

e 


dN t . 


-{a i +p i x+ r t y) 


2 A' 


\Pi] 


(3.111) 

(3.112) 

(3.113) 


Los valores de a i , /L y y k son constantes una vez que se conoce las coordenadas de los 
nodos. 

<X> ; , ® y O k son independientes de las coordenadas espaciales, por lo que su derivada 
tiene un valor constante. 

En gradiente constante en cualquier elemento significa que debe usarse un gran numero de 
elementos para aproximar un rapido cambio de (j> . 

EJEMPLO 3.7 

Evaluar las funciones de forma y calcular el valor de la presion en el punto A. Los valores 
nodales son d> ; = 40 kg/cm 2 ; ® . = 34 kg/cm 2 y ® t = 46 kg/cm 2 . El punto A esta situado en 

(2,1.5). 



Nodo 

X 

Y 

i 

0 

0 

./ 

4 

0.5 

k 

2 

5 


Fig. 3.22 Valores y coordenadas nodales del ejemplo 3.7. 


a t = x } y k ~x k yj = 4(5)-2(0.5) = 19 
a j = Vk, -x,y k = 2(0)-0(5) = 0 
a k = x t y j - x j y i = 0 (0.5) - 4 (0) = 0 
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P i =y ] -y k = 0.5-5 = - 4.5 
P j =y k -y i = 5-o = 5 
P k =y-yj = 0-0.5 = -o .5 


U = ** -x y =2-4 = -2 


y j =x i -x k =0-2 = -2 
r k =x j -x i = 4-0 = 4 

Determinante: 



1 x t y, 


1 0 0 

2 A e = 

1 *7 ^ 

= 

1 4 0.5 


1 x t y k 


1 2 5 


20-1 = 19 


2aT =19 


Sustituyendo en las funciones de forma: 

w ‘ =ttt[-° 5x+4v ] 


Se debe comprobar que Af + Ay + A^ = 1 


— [l9-4.5x-2y + 5x-2_y-0.5x + 4y] = 1 


La ecuacion para interpolacion es: 


^ = [l9-4.5x-2y]^O, + [5x-2y]^O /+ [-0.5x + 4y]^a) A 


en.4 (2, 1.5) 


k i \ 


v 19y 


r 5 


*= T7T K 40 ) + T7T K 34 ) + T7T K 46 ) 


19 


v 19. 


(f> = 39.40 kg/cm 2 


Se satisface que las funciones de forma varian linealmente a lo largo de los lados entre su 
nodo y los otros dos nodos, es decir, A) varia linealmente a lo largo de los lados ij y ik. Una 
funcion de fonna es cero a lo largo del lado opuesto a su nodo, esto es Ni es cero a lo largo 
del lado jk. 
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EJEMPLO 3.8 

Determinar la localizacion de la llnea de contomo de 42 kg/cm 2 para un elemento triangular 
del ejemplo anterior. 


<t>; = 40 



<*>£ = 


46 


/ = ( 0 , 0 ) 



Fig. 3.23 Elemento triangular del ejemplo 3.8. 
Se puede ver que la presion para 42 kg/cm 2 intercepta los lados ik y jk. 
Para el lado jk 


46-42 _ 2-x 
46-34 ~ 2-4 

4 _2-x 
12 ~ -2 


x = 2.67 cm 


46-42 _ 5-y 
46-34 ” 5-0.5 

y = 3.5 cm 

en el lado jk: x = 2.67 cm 

y = 3.5 cm 

Para el lado ik: 

40 - 42 _ 0 - a 
40-46 “ 0-2 


x = 0.666 cm 


40-42 _ 0 -y 
40-46 “ 0-5 


y= 1 .667 cm 


3.5 ELEMENTO RECTANGULAR BILINEAL 


El elemento rectangular bilineal tiene una longitud 2b y una altura 2a. Los nodos son i,j, k 
y m. El nodo i siempre se ubica en la esquina inferior izquierda. 
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<&n 


* 

m 


k 

> — 

i 


j 


2b 



t 

r 


'm 

k 



L 

J 

i 

» 

j 

4 

s 

p — 


Fig. 3.24 Elemento rectangular bilineal. 

La ecuacion de interpolacion en coordenadas locales s y t es: 

<j) = Cj + c 2 s + c 3 t + c 4 st (3.114) 

esta funcion es lineal en s a lo largo de cualquier linea constante t y lineal en t a lo largo de 
cualquier linea constante s. 

El sistema coordenado local s-t tiene su origen en el nodo i. Evaluando la funcion de 
interpolacion en los nodos tenemos: 



o 

ll 

o 

II 

<< 

0,- = c i 

(3.115) 


s-2b\ t - 0 ; 

0 j = Cj + 2 bc 2 

(3.116) 


s = 2b ; t = 2a ; 

0, = Cj + 2 bc 2 + 2 ac 3 + 4abc 4 

(3.117) 


cs 

CN 

II 

o 

II 

0»=c 1 +2ac 3 

(3.118) 

De donde: 

Cl = O, 


(3.119) 

En 

O J = 0, + 2 bc 2 


(3.120) 

Despejando 

0-0, 

= 

2b 


(3.121) 

En 

=®, + 2ac 3 


(3.122) 

Despejando 

0 -0 

c, = 

2 a 


(3.123) 


Sustituyendo c/, C 2 y c? en: 


0, =d>, +2 b 


2b 


+ 2 a 


O -O. 

m i 


2 a 


+ 4 abc. 


(3.124) 
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Despejando: 


— [~0 -O . + 0, - 0 1 

L 1 j k 


4 ab 


Sustituyendo en 

(j) = Cj + c 2 s + c 3 t + c 4 st 


<f> = 0, + 


0 ,- 0 / 


2b 


s + 


Tomando terminos comunes: 

st t 


</> = ■ 


1 - — 
2b 


+ 


4 ab 2 a 
Que se puede escribir: 


> 0 , +< 


0 -0 

m i 

2 a 


st 


(3.125) 

(3.126) 


0 127 ) 

4 ab 


2b 


4 ab 


0 , + 


st 

4 ab 


> 0 A +< 


2 a 


st 


> 0 ., 




, s 

r t 

0,+ 

s 

r 1 1 

0, + 

st 

0, + 

’ t ’ 

, s 

1 — 

i — 


i — 

1 

2b J 

2 a _ 


\_2b\ 

2 a 

J 

4 ab _ 


2a _ 

2b\ 


Donde: 


0 


^ = JV I .0,.+iV.0.+iV i 0,+7V m 0„ 


4 ab 

(3.128) 

(3.129) 

(3.130) 


N,= 

N J = 

N k = 
N m = 




s 

2b 



st 

4 ab 


t 

2a 



(3.131) 

(3.132) 

(3.133) 

(3.134) 


Las funciones de forma de un elemento rectangular bilineal tienen propiedades similares a 
las de un elemento triangulo: 

1.- Cada funcion de forma varia linealmente a lo largo de sus ejes entre su nodo y dos 
nodos adyacentes, por ejemplo, N, varia linealmente a lo largo de los nodos ij y im. 

2 - Cada funcion de forma es cero a lo largo de los lados cuyo nodo no toca, por ejemplo, N, 
es cero a lo largo de los lados jk y km. 


3.- La variacion lineal de <j> a lo largo de los ejes del elemento rectangular y un lado de un 
elemento triangular significa que esos dos elementos son compatibles y pueden ser usados 
adyacentemente uno de otro. 



Fig. 3.25 Compatibilidad de elementos. 


150 


FUNDAMENTOS BASICOS DEL ELEMENTO F1NITO Hernandez Barrios Hugo FESA-UNAM 


Otro sistema coordenado popular en el que tiene su origen en el centro del elemento. 


b b 



Fig. 3.26 Sistema coordenado con origen en el centro del elemento. 

La relacion entre ambos sistemas coord enados es: 

s = b + q 
t = a + r 

Las funciones de forma son: 


' 4 

"l-«" 

b_ 

\- r - 
a _ 

(3.135) 

Nj = - 

1 4 

1 l 

+ 

1 1 

\- r - 

a 

(3.136) 

N k =- 

k 4 

1 1 ^ 

+ 

J 1 

\ + r - 
a _ 

(3.137) 

N ", ~ ^ 

i 

1 

\ + V - 

a 

(3.138) 


estas funciones de forma son utiles debido a que permiten que un rectangulo se deforme en 
un cuadrilatero general. 

Para el sistema coordenado natural, la relacion entre coordenada es: 


Donde qy r son coordenadas locales. 


'=! 

(3.139) 

r 

v = - 

(3.140) 


a 
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Fig. 3.27 Coordenadas locales q y r. 


Este sistema coordenado tiene su origen en el centro del elemento y es conveniente en el 
caso de realizar integracion numerical 




n = -i 

: = -i ; = -i 

Fig. 3.28 Sistema coordenado con origen en el centro del elemento. 



- 1<^<1 
-1 < 77 < 1 



(3.141) 


(3.142) 


(3.143) 


(3.144) 


Una linea de contomo en un elemento rectangular es generalmente curva, la interseccion de 
la linea de contomo con los ejes se puede obtener por medio de interpolacion lineal, el 
tercer punto se pude obtener haciendo sot igual a cero en la ecuacion de la funcion de 
forma y resolver para otro valor coordenado. 

EJEMPLO 3.9 

Detenninar los tres puntos de la linea de contomo para 50° C del elemento rectangular 
mostrado, los valores nodales son: 

cp . = 42°C = 56°C 

Oj = 54°C O m = 4 6°C 
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Fig. 3.29 Elemento rectangular del ejemplo 3.9. 


Longitudes de los lados: 

26 = Xj - x t 

26 = 8-5 = 3 ; 26 = 3 

2 a=y m ~y i 

2a = 5-3 = 2 ; 2a = 2 

Sustituyendo en las funciones de fonna (3.13 1,3-132, 3. 133, 3.134): 



La linea de 50°C intersecta los nodos ij y mk. Por lo que suponemos los valores de t y 
calculamos los de s. 

1) A lo largo de los nodos ij ;t = 0 



N k = 0 
N m = 0 


d>= 1-- 0.+ - O = 1-- (42)+ - (54) = 50 

L 3J ' [_ 3 J 1 L 3J L3J 

42 54 

42 s + — s = 50 

3 3 
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45 = 8 


2 ) 


A lo largo del lado mk , t = 2a 



N k = 
N = 

m 



St _ 5(2) _ 5 


6 6 3 


~ t~ 

"l-^' 



' 2 

\- S - 

_2_ 

3_ 


_2_ 

3_ 


2 



5 = 2 
t = 0 


<P 





56 

T 


5 - 


— 5 + 46 = 50 
3 


(56) 


+ 


1 - 


(46) = 50 


5 = 1.20 
t = 2.0 


Para el tercer punto, suponemos que t = a = 1 (centro del elemento) 


0 = 
50 = 


1--1 
3 J 

'i-ll 
2 _ 

®,- + 

5 

_3_ 

1- 

"l-^" 

3_ 

r r 
L2_ 

(42) + 

5 

_3_ 

" 1 " 

_2_ 


o ; + -o A + 

6 


1 -* 


4), 


(54) + — (56) + 
o 


1 -^ 


(46) 


50 = 21- — 5 + 95 +— 5 + 23-— 5 
3 6 3 

50 = 44 + — 5 
3 


5 = 1.64 
t = 1 



(2,0)d(7,3) 


Fig. 3.30 Linea de contorno del ejemplo 3.9. 
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3.6 DEFORMACIONES PLANAS 

El estado de deformaciones planas ocurre en miembros que no son libres de expandirse en 
la direccion perpendicular al piano de aplicacion de la carga. Si suponemos que la carga 
aplicada se encuentra en el piano x-y, entonces los desplazamientos w que se presentan en 
la direccion z son cero, y los desplazamientos u y v son funcion unicamente de x y y. 

Este conjunto de desplazamientos hace que las deformaciones y x _ = y = 0 y que s z sea 
constante o cero. 

El tensor de deformaciones es: 

N 


e * 2 Yxy ° 

\r*y £ y 0 

0 0 0 


Son ejemplos de defomiaciones planas: un tunel, una tuberia a presion interna, una presa de 
gravedad, etc. 


3.7 ESFUERZOS PLANOS 


Otro caso de interes se presenta si restringimos todos los esfuerzos actuando en un solo 
piano. Lo anterior es similar al caso de deformaciones planas, aunque los esfuerzos a x son 

cero, los valores de s_ no son necesariamente constantes. Por tanto, las secciones 

transversales pueden no permanecer planas y el estado de esfuerzos pianos se convierte en 
un problema de defomiaciones tridimensionales. 

Un estado de esfuerzos pianos existe cuando un cuerpo elastico es delgado y no hay cargas 
aplicadas en la direccion paralela al espesor. Cuando la carga es aplicada en el piano x-y las 
componentes de los esfuerzos asociadas con las direcciones del espesor a_, r xz y r yz son 

muy pequenas y por tanto se pueden suponer igual a cero. 

Son ejemplo de esfuerzos pianos: una viga peraltada, una pieza sometida a tension, etc. 


El tensor de esfuerzos es: 


H 



xy 

0 

0 (7 _ 


3.8 RELACIONES ESFUERZO-DEFORMACION 

Para un material isotropico, la matriz constitutiva que relaciona los esfuerzos con las 
defomiaciones, para el caso de Esfuerzos pianos es: 
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cr 


(1 + *')(!-•') 

y para el caso de deformaciones planas: 


1 v 
v 


0 0 


xy 


(l + v)(l-2v) 


0 

1 0 
1 — v 


1 — V 
V 

0 


/r 


V 0 

l-v 0 
l-2v 


0 


L Yxy 


Ambas expresiones se denotan: 

donde [D] es la matriz constitutiva del material y siempre contiene tenninos constantes 
que son reflejo de las propiedades del material. 


3.9 ELEMENTO TRIANGULAR 




Fig. 3.31 Componentes de desplazamiento. 



El desplazamiento de un punto A, para el caso de esfuerzos pianos, se puede representar por 
los componentes de sus desplazamientos en las direcciones x-y ; u y v ; asi cada nodo tiene 
dos componentes de desplazamiento y el numero total de componentes de desplazamiento 
es de 6, lo anterior se puede observar en la figura (3.31). 

Como u y v son independientes uno de otro, los coeficientes de las funciones de 
interpolacion no deben ser los mismos. Asi podemos proponer un polinomio de 
interpolacion lineal para cada componente de desplazamientos. 


Se habia visto que (j) = a, + a 2 x + a 3 y ; donde (j) es el valor nodal o incognita. 


por nodo: u = + a 2 x + a 3 y 

v = a 4 + a 5 x + a 6 y 


(3.145) 

(3.146) 
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Los desplazamientos son independientes, se pueden expresar: 



a j 


l x y 0 0 
0 0 0 1 a 


0 

y 


a 2 

a 3 

a 4 


a. 


l«6 


Para el total de nodos: 


v 


‘1 

L 

Li 

0 

0 

0 " 


a x 

L 


0 

0 

0 

i 

L 

y t 


a 2 

U] 


1 

X J 

yj 

0 

0 

0 

< 

a 3 

> 

v j 


0 

0 

0 

1 

X J 

y j 


a 4 

u k 


1 

x k 

y k 

0 

0 

0 


a 5 

v k . 


0 

0 

0 

1 

** 



a 6 


( 3 . 147 ) 


( 3 . 148 ) 


Lo anterior implica que los valores de a,, a 2 ...a e son los mismos para cada nodo y por 
tanto el elemento es indeformable en sus lados que temporalmente podemos escribir: 


' ■> 

u < 


“1 

*,■ 

Li 

0 

0 

0 " 

u j 


1 

X J 

yj 

0 

0 

0 

u k 


1 

*k 

y k 

0 

0 

0 

L 


0 

0 

0 

1 

X, 

Li 

v j 


0 

0 

0 

1 

Xj 

y j 

v k. 


0 

0 

0 

1 

*k 

yk. 


Invirtiendo la matriz y despejando: 

{ a 1 = MM 
M = M 1 {«"} 

A A 

coef. incognitas grados de libertad 


( 3 . 149 ) 


( 3 . 150 ) 

( 3 . 151 ) 

( 3 . 152 ) 
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(*3 

- x k yj ) 

(• x k y t ~ x ,y k ) 

< V,V: -Xj-yJ 

0 

0 

0 

«/ 


( yj 

-y k ) 

(. y k -y , ) 

(y t-yj) 

0 

0 

0 

u i 

1 

( x k 

-Xj) 

( x i~ x k) 

(Xj-Xi) 

0 

0 

0 

u k 

2 A 


0 

0 

0 

C ~x t yj) 

(• x k y t ~ x iy k ) 

I-'',.'',. -Xjy t ) 

y 



0 

0 

0 

( y } -y k ) 

tyt-yi) 

(y i-yj) 

v , 



0 

0 

0 

(x k -Xj) 

(x, -x k ) 

(Xj-X,) 

3. 


Donde 


A = 


x, 

x 


1 X, 


y t 

yj 

y k 


area del triangulo i-j-k 


( 3 . 153 ) 

( 3 . 154 ) 


Organizando: 


0\ 


{xjy.-x^j) 

0 


0 

(xyj-xjy] 

1 0 

U i 

a. 


iyj-yt) 

/ \ 

0 

04-3/) 

0 

iy-yj) 

/ \ 

0 

y 

a, 

1 


0 


0 

\ x j- x i) 

0 

u j 

a, 


0 

(.Xjy k -x k yj) 

0 

{x k y-Xiy k ) 

0 

(vi-Vl) 

v j 

a 5 


0 

( yj-Vk ) 

/ \ 

0 

(y k -y,-) 

0 

(>;->•,) 

/ \ 

u k 

_«6 J 


0 


0 

( x i-x k ) 

0 

(3 -3) _ 

3. 


( 3 . 155 ) 


Como: 


{«} = [c e ] { a e } 
u = a { + a 2 x + a 3 y 
v = a 4 + a 5 x + a 6 y 
du 


dx 

dv 

dy 


■ = a. 


= a. 


du dv 
dy dx 


+ — = a 3 + a 5 


( 3 . 156 ) 

( 3 . 145 ) 

( 3 . 146 ) 

( 3 . 157 ) 

( 3 . 158 ) 

( 3 . 159 ) 


Calculo de deformaciones, caso esfuerzos pianos 
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M = 


e , 

s 

r„ 


du 


a x 

dx 


“0 1 0 0 0 0“ 

a 2 

dv 



a , 


= 

0 0 0 0 0 1 

j 

~dy 



a 4 



0 0 10 10 


du dv 


L - 1 

a 3 

1 


D 

dy dx 


_a 6 _ 


(3.160) 


M =[(?]{«} 

M=[0]M'V} 




(y,-y k ) 


o 


o 


iyk-yd 0 (y,-yj) 

0 (x k -Xj) 0 (*,--**) 0 (Xj-x,) 

(X k -x.) (y,-y k ) (x i -x k ) {y k -y t ) (x.-x t ) (y,-*,) 


W =[«](“1 

Otra forma 6 bien como: u = N l u l + N 2 u 2 + N 3 u 3 

v = 7V 1 v 1 +7V 2 v 2 +7V 3 v 3 


Por nodo 


(3.161) 

(3.162) 

a e J (3.163) 

(3.164) 

(3.165) 


H = - 


M 


0 

n 2 

0 

n 3 

0 “ 


0 

N, 

0 

n 2 

0 

7V 3 _ 


L v 3 


(3.166) 


Con TV. = — — (a, + bx + c,y ) i = 1,2,3 

i v 1 1 l,/ y 

Como anteriormente se definio 

n ‘ = Ya 7 + ^‘ x + r,} ^ i = l ’ 2 ’ 3 


(3.167) 


(3.168) 


du 

dN , 

dN , 

dN 3 

= 

= M, 

+ — - 

M, H i 

dx 

5x 1 

dx 

2 cbc 

dv 

57V, 

dN , 

57V , 

= 

= — L v, 

+ 

-V, H 1 

dy 

dy 

5y 

5y 


(3.169) 

(3.170) 


159 



FUNDAMENTOS BASICOS DEL ELEMENTO FIN1TO Hernandez Barrios Hugo FESA-UNAM 


du dv dN. 
x = — h = — -u, + ■ 

/ XV ^ ^ /-.l 


dN, 


dy dx dy 


dx 


-Vi+- 


dN , 


dy 


-U-, + ■ 


dN 


dx 


-v 2 +- 


dN 


dy 


- U 3 + ■ 


dN 


dx 


du 


dN, 

0 

dN , 

0 

dN 


u \ 




3 

0 


dx 


dx 


dx 


dx 


v i 

dv 


0 

dN, 

0 

dN , 

0 

dN 

U 2 


= 

2 

3 





dy 


dy 


<3y 

V 2 

du dv 
1 


dN, 

dN, 

dN 2 

dN 2 

cW 3 


u 3 

dy dx 


_dy 

dx 

dy 

dx 


dx 

_v 3 _ 


donde: 


M =[*]{«*} 

[b} = [b v b 2 ,b,] 


B, = 


dN, 

dx 


0 


dN, 


dy 

dN, dN, 


dx dy 


Matriz de deformaciones del nodo i. 


N i=-^[ a i +b i x + c iy] 


N 2=Y~ l l a 2 + b 2 X + C 2}’\ 


N 3 = ~~[ a 3 + b 3 x + c 3 y] 

A xi. 


[*] = 


2 A e 


b x 0 b 2 0 b 3 0 

0 Cj 0 c 2 0 c 3 

c, b ] c 2 b 2 c 3 b 3 


y por consiguiente: 


B = 


2 A e 


\ 0 

0 c ; 
c, b, 


i — 1,2,3 


(3.171) 

(3.172) 

(3.173) 

(3.174) 

(3.175) 

(3.176) 

(3.177) 

(3.178) 

(3.179) 

(3.180) 
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[#] es una matriz de constantes y por tanto las deforaiaciones son constantes a traves del 
elemento. 


b x = ( yj-y k ) G = (** ~Xj) ( 3 . 181 ) 

b 2 = ( y k -yt ) C 2 = (x, - x k ) ( 3 . 182 ) 

63 = O',- - y } ) c 3 = (Xj -x t ) ( 3 . 183 ) 

Campo de esfuerzos: 

[<r] = [ D M = [ D ][ B ]{°'} ( 3 . 184 ) 


Los esfuerzos tambien son constantes a traves del elemento. Si existen defonnaciones 
iniciales y/o esfuerzos iniciales. 


Para: 


donde: 


H = [D][5]{^}-[D][ Go ] + K] (3-185) 





0 

d' d^2 0 
0 0 d 33 


(3.186) 


(3.187) 


Esfuerzos pianos 


d 

d 

d 

d 

d 


11 


12 


21 


22 


33 


E 

(1-v 2 ) 

vd n 

vd n 

d] 1 

E 

2(1 + v) 


G 


Defonnaciones planas 
E(l-v) 


d\ 1 — 


(1 + v )(l - 2v) 


d ^ d 1 2 


t/p d-y^ d |2 


1-v 




2(1 + v) 


= G 


(3.188) 
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La matriz elemental de rigidez es: 


[*■] = /, [*]>][*]*■ 

*' v 01 


Para un elemento en particular [#] es constante y para un espesor t, constante en el 
elemento. 


[k e ] = [B] T [D][B]tjdxdy 
[k-} = At[B] T [D\[B] 


A es el area del elemento triangulo. 


\[b\[D]{B]Jv = A,[B]{D]IB] 



t 

4A 


b, 0 

0 Cj 

b 2 0 

0 c n 


b 3 0 

0 c 3 


d\\ 

d\2 

0 " 

h \ 

0 

b 2 

0 

b 3 

0" 


d 2l 

d 22 

0 

0 

c 

0 

c 2 

0 

c 3 

(3.189) 

L 0 

0 

d 33 



c 2 

b 2 

C 3 

b 3 













EJEMPLO 3.10 


Determine los desplazamientos, deformaciones y esfiierzos para la viga en cantiliver 
cargada como se muestra en la figura (3.32). 


-r 



E = 22135947 tm 2 
v = 0.25 
t = 0.1m 


Fig. 3.32 Viga en cantiliver. 


Se divide la viga en dos elementos triangulos con el numero de nodos siguientes: 


1 3 



4 


Fig. 3.33 Topologia de la viga. 


Coord enadas de los nodos: 
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Nodo 

X 

y 

1 

0 

0.35 

2 

0 

0 

3 

6 

0.35 

4 

6 

0 


ELEMENTO 

Nodo 

i 

i 

k 

1 

1 

2 

4 

2 

1 

4 

3 


Condiciones de frontera: 


M, = v, = 0 
u 2 = v 2 = 0 
F = -2 

r 3y ^ 

Incognitas: 

Desplazamientos « 3 ,v 3 ,w 4 ,v 4 

Reacciones R u ,R tv ,R 2y ,R 2x 


Area de los elementos 




"1 

0 

0.35" 


A. 1 - 

1 


y i 

elemento 1 

2A 

1 

0 

0 

= 2.1 

1 

x 2 

V 2 







2 





1 

6 

0 



1 

X 4 

t 4 


elemento 2 

Matrices elementales: 

Elemento 1 


4 =1.05 
A 2 =1.05 

i j 

1 2 

(0,0-35) (0,0) 


k 

4 nodos 

(6,0) coordenadas 


b f = Vj -y k = 0 
h j =y k -y t =-0.35 
b k =y t - y, =0.35 


c i =x k ~x J =6 
Cj=x~x k = -6 
c, = x, — X- =0 

/C J l 
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[B] 


0 0 -0.35 
0 6 0 
6 0-6 


0 0.35 0 

-6 0 0 
-0.35 0 0.35 


Esfuerzos pianos 


E = 2213594 (T/m 2 ) 


dl 1 = E / (1 - v 2 ) = 2361166.93 

dl2 = vdll = 590291.73 

d21 = vdll = 590291.73 

d22 = dll = 2361166.93 

d33 = E / 2 (1 + v) = 885437.60 


Id] 


d 11 d 12 0 

dll dll 0 
0 0 d33 


Defonnaciones planas 


Id] 


2361166.93 590291.73 0 

590291.73 2361166.93 0 

0 0 885437.60 


[k e ] = [B] T [D][B]tA 


k c \ = (0.1)(1.05) 


0 

0 

2.86 

0 

2.86 

0 

-0.17 

0 

-2.86 

0 

-2.86 

-0.17 

0.17 

0 

0 

0 

0 

0.17 


2361166.93 590291.733 
590291.733 236116693 
0 0 


0 

0 

885437.6 


0 

0 

2.86 


0 -0.17 0 0.17 0 

2.86 0 - 2.86 0 0 

0 -2.86 -0.17 0 0.17 


U | 

v i 

u 2 

V 2 

u 4 

V 4 


76046551 

0 

-76046551 

-4520250 

0 

4520250" 

U | 

0 

202790651 

-3013500 

-202790651 

3013500 

0 

v i 

-76046551 

-3013500 

76763017 

7533749 

-716496 

-4520250 

U-, 

-4520250 

-2027906 S 1 

7533749 

20305917 

-3013500 

-268686 

V 2 

0 

3013500 

-716496 

-3013500 

716496 

0 

u 4 

4520250 

0 

-4520250 

-268686 

0 

268686 

V 4 
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Elemento 2 


j k 


1 4 3 

(0,0.35) (6,0) (6,0.35) 


nodos 

coordenadas 


b, = -0.35 
bj= 0 
b k =0.35 


c,.=0 

= - 

c k = 6 


[B] 


1 

2A 


0.35 

0 

0 


0 0 0 0.35 0 

0 0-60 6 
-0.35 -6 0 6 0.35 


[rl=[B] T [D][B]tA 


ll | V j M 


V 4 U 3 V 3 



7164.96 0 

0 2686.86 

0 4520250 

30135.00 0 

- 7164.96 - 45202.50 

- 30135.00 - 2686.86 


0 30135.00 

45202.50 0 

760465.51 0 

0 2027907.17 

- 760465.51 30135.00 

- 45202.50 - 2027907.17 


- 7164.96 

-4520250 

- 760465.51 

- 30135.00 

767630.47 

75337.49 


- 30135.00 

- 2686.86 

- 45202.50 

- 2027907.17 

75337.49 

2020594.03 


U , 

V 1 

V 4 

m 3 

V 3 


Matriz de rigidez total: 


M 


U J 

v l 

U 2 

V 2 

' 76763Q47 

0 

-760465.51 

-4520250 

0 

203059317 

30135.00 

-202790631 

-760465.51 

-30135.00 

76763Q47 

75337.49 

-4520250 

-202790631 

75337.49 

203059317 

-7164.96 

-4520250 

0 

0 

30135.00 

-2686.86 

0 

0 

0 

75337.49 

-7164.96 

-30135.00 

75337.49 

0 

-4520250 

-2686.86 


U 3 

v 3 

U 4 

V 4 


-7164.96 

-30135.00 

0 

75337.49 

«i 

-4520250 

-2686.86 

75337.49 

0 


0 

0 

-7164.96 

-4520250 

« 2 

0 

0 

-30135.00 

-2686.86 

V 2 

76763Q47 

75337.49 

-760465.51 

-30135.00 

“3 

75337.49 

203059403 

-4520250 

-2027907.17 

V 3 

-760465.51 

-4520250 

76763Q47 

0 

« 4 

-30135.00 

-2027907.17 

0 

203059403 

'’4 


M| = 0 

Vj = 0 
= 0 
v 2 = 0 


Relacion fuerza desplazamiento [^]{«} = {i 3 } 


165 


FUNDAMENTOS BASICOS DEL ELEMENTO FIN1TO Hernandez Barrios Hugo FESA-UNAM 


" 76763047 

7533749 

-76046551 

-3013500 ' 

u 3 


" 0 " 

7533749 

203059403 

-4520250 

-202790717 

V 3 


-2 

-76046551 

-4520250 

76763047 

0 

U 4 


0 

3013500 

-202790717 

0 

203059403 

y 4 _ 


0 


Resolviendo: 


u 3 


"- 0.0000727 

-0.0000124 

-0.0000713 

-0.0000113" 

' 0 ' 


" 0.0000248 ' 

V 3 


-0.0000124 

0.0003713 

0.0000095 

0.0003706 

-2 


- 0.0007426 

u 4 


0.0000713 

0.0000095 

0.0000725 

0.0000106 

0 


-0.0000191 

y 4 _ 


-0.0000113 

0.0003706 

0.0000106 

0.0003704 

0 


-0.0007412 


0.0000248 



1 


0.0000191 


0.0007426 


0.0007412 


Fig. 3.34 Desplazamientos de la viga en cantiliver. 


Desplazamiento Vertical Promedio 


0.0007426 + 0.0007412 
2 


0.0007419 


S _PL 3 _ (2)(6) 3 (12) 

3 El 3(221 3594)(0. 1)(0.35) 3 

0.18 + 0.0007419 


0.18 


Que puede mejorarse la aproximacion refinando la malla. 
Esfuerzos: [cr] = [Z>][5]{a e | 
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Elemento 1 


07 


"2361166933 

5902917333 

0 


= 

5902917333 

2361166933 

0 

.V 


0 

0 

885437.6 


o 







0 



0 

-0.17 

0 

0.17 

0 ' 

0 


' -7.6667 ' 

2.86 

0 

-2.86 

0 

0 

0 

— 

-1.9166 

0 

-2.86 

-0.17 

0 

0.17 



-111.5687 






-0.0000191 




-0.0007412 


Elemento 2 













0 

0 



<7 X 


‘2361166933 

5902917333 

0 

-0.17 

0 

0 

0 

0.17 

0 ' 



‘ 7.5911 " 












-0.0000191 



a 

= 

5902917333 

2361166933 

0 

0 

0 

0 

-2.86 

0 

2.86 

-0.0007412 

= 

-6.9654 

(7 


0 

0 

8854376 

0 

-0.17 

-2.86 

0 

2.86 

0.17 



-0.6091 

L xyj 











0.0000248 


- 












-0.0007426 





Fig. 3.35 Esfiierzos de 



Ahora se mostrara la manera en la cual refinando la malla se mejora la aproximacion de los 
resultados del ejemplo de la viga en cantiliver. 

Esto se hace en el Sap con el elemento Shell (cascara), el cual se usa para modelar cascaras, 
membranas y comportamiento de placas en estructuras planas y tridimensionales. Este 
elemento es una formulacion de tres a cuatro nodos que combinan el comportamiento 
membrana y flexion. 


VIGA EN CANTILIVER DEL EJEMPLO 3.10 [P.154.] 
ELEMENTO TRIANGULAR 

Primero veamos los resultados con la membrana de 4 nodos: 
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X ■ ■ .1 w 

File Edit View Define Bridge Draw Select Assign Analyze Disglay Design Options Help 

1 D & | B ^ I --T r. | f | [a >1 fT# @ jS 0 & | H | 3-d n. a k '» CT to | O- ■& | H j . | 


yi*js 




nhH- . 


£ t- H IE. * 


it J 


-H 

: a 

l< 


X Deformed Shape (DEAD) 


□EMI X X-Z Plane 


□E3B; 



Right Click on any joint for displacement values 


Start Animation 1^;^! GLOBAL j|Ton.m.C * 


Como podemos observar, los resultados de los desplazamientos son los mismos o muy 
parecidos a los realizados manualmente. 
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Manualmente [m] 

u 3 = 0.0000248 
v 3 = -0.0007426 
m 4 =-0.0000191 
v 4 =-0.0007412 


En el Sap [m] 

u 3 = 0.0000263 
v 3 =-0.0008154 
u 4 = -0.0000205 
v 4 =-0.0008148 


Ahora veamos los resultados de los esfuerzos 



Manualmente [T] 


En el Sap [T] 


Elemento 1 

Elemento 2 

Elemento 1 

Elemento 2 

(j x = -7.6667 

a x =7.5911 

cr x = —12.2259 

a x = 12.6363 

<t v =-1.9166 

ct v = -6.9654 

<7 Y =-3.0564 

<7 y =-3.9412 

cr =-111.5687 

xy 

=-0.6091 

<7 =-93.7723 

X)/ 

a ^ = -29.5200 
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MIN=-1 2.226, MAX=1 2.636, Right Click on any Shell Element for detailed diagram 


O | ={> (GLOBAL ^]|Ton. m.C - 


A continuation veamos los resultados con una membrana de 10 nodos: 



Veamos como los resultados se empiezan a aproximar 
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Desplazamientos [m] 

« 3 = 0.00033 
v 3 = -0.00796 
u 4 =-0.00031 
v 4 = -0.00796 


SAP2000 vlO.0.1 Advanced - membrana3 




File Edit View Define Bridge Draw Select Assign Analyze Display Design Options Help 

|| D Bg| H S> | « r, f |a I ► @ @ @@ & HU 3<1 W kz yz n» Cf iftf | O’ Q \ ?u 13 | 

II » t 

ff 



MIN=-305.476, MAX=307.787, Right Click on any Shell Element for detailed diagram 


O | ^ | GLOBAL ]£j|Ton.m,C 


Esfuerzos [T] 


Elemental 

Elemento 8 

o x = -305.4757 

a x = 39.7258 

<j v =-76.3689 

u y =-19.3566 

ct =-311.0863 

xy 

= -26.2062 

xy 


A continuation veamos los resultados con una membrana de 26 nodos: 
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$£ ;AP2000 vIO.O. 1 Advanced Tiemlyanaf 

File Edit View Define Bridge Draw Select Assign Analyze Display Design Options Help 


iHQ 


a . n "n H - . 



Right Click on any joint for displacement values 


Start Animation 


~\<P | => (GLOBAL -j|Ton,m,C - 


Observemos que los resultados se siguen aproximando 


Desplazamientos [m] 

u 3 = 0.0015 
v 3 = -0.0352 
u 4 = -0.0015 
v 4 = -0.0352 


Ahora veamos los resultados con una membrana de 38 nodos: 
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Desplazamientos [m] 

« 3 = 0.003 
v 3 = -0.081 
u A = -0.003 
v 4 =-0.081 

Como se ha demostrado, el elemento Shell es muy efectivo en el analisis estructural y en 
este ejemplo se ha utilizado la fomia triangular, sin embargo la fonnulacion cuadrilatera es 
la mas exacta de las dos[l]. 


[1] Manual Sap2000 con ejercicios en espanol, P. [55] 
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CONCLUSIONES 


Historicamente las primeras aplicaciones del metodo del elemento finito fueron el analisis 
de defonnaciones y esfuerzos en cuerpos y estructuras grandes, ademas de problemas de 
flujo de calor, campos magneticos, etc. Actualmente existen muchos programas 
comerciales que utilizan como base este metodo. 

Conclusiones Particulares 

Existen diversos metodos para resolver ecuaciones diferenciales, dentro de los cuales se 
encuentra el metodo del elemento finito, que basado en una formulacion integral (metodo 
de residuos pesados), se utiliza el metodo de Galerkin como mejor opcion ya que es mas 
facil de programar porque se puede colocar la funcion de peso igual al polinomio de 
aproximacion. 

Los sistemas coordenados se clasifican en: sistemas locales, sistemas naturales y sistemas 
de area. Estos sistemas nos facilitan la solucion de las integrales que necesitan ser resueltas 
en el metodo del elemento finito. 

Dentro de los elementos unidimensionales y bidimensionales existen una gran cantidad de 
elementos. Estos elementos pueden ser resueltos por muchos metodos, los cuales derivan 
las ecuaciones diferenciales para obtener la solucion. El metodo del elemento finito utiliza 
la integracion para resolver los elementos y a partir de encontrar la matriz de rigideces se 
encuentra la solucion de la misma manera. 

Existen varias ventajas al utilizar este metodo: 

1. - La solucion numerica de ecuaciones diferenciales se puede formular en terminos de una 
integral. 

2. - El metodo utiliza fimciones continuas para la aproximacion de una cantidad o 
cantidades desconocidas. 

3. - Como caso particular, al utilizar el metodo de Galerkin se puede colocar la funcion de 
peso igual al polinomio de aproximacion con lo cual es mas facil de programar. 

4. - Una vez encontrada la matriz de rigideces, se es libre en el metodo a elegir para la 
solucion de la misma. 

Ademas de las ventajas del metodo se pueden mencionar las aportaciones de este trabajo: 

1. - Se hizo una recopilacion de los conceptos basicos del elemento finito, desde la 
formulacion integral (metodo de residuos pesados) hasta la integracion numerica. 

2. - Se demostro que de un tema complicado como lo es el elemento finito, se pueden 
plantear y solucionar problemas de ingenieria civil con la misma facilidad que lo hacen 
otros metodos. 
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Conclusiones Generales 

Los ejemplos desarrollados manualmente corresponden a los del programa Sap 2000 V.10, 
utilizado en este trabajo, con lo que de alguna manera se puede llegar a demostrar que este 
tipo de programas comerciales realizan el mismo procedimiento analltico. 

En este trabajo se mostraron las bases de algunos de los elementos como el elemento barra, 
elemento viga, elemento triangular, etc. Dejando a futuro otros elementos por desarrollar. 
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